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Abstract 

We give an introduction to three topics in lattice gauge theory: 

I. The Schrddinger functional and 0(a) improvement. 

0(a) improvement has been reviewed several times. Here we focus on explaining the 
basic ideas in detail and then proceed directly to an overview of the literature and our 
personal assessment of what has been achieved and what is missing. 

II. The computation of the running coupling, running quark masses and the extraction 
of the renormalization group invariants. 

We focus on the basic strategy and on the large effort that has been invested in under- 
standing the continuum limit. We point out what remains to be done. 

III. Non-perturbative Heavy Quark Effective Theory. 

Since the literature on this subject is still rather sparse, we go beyond the basic ideas 
and discuss in some detail how the theory works in principle and in practice. 
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1 Introduction 



Lattice QCD, the subject of this school, is our prime tool to make quantitative pre- 
dictions in the low energy sector of QCD. Also connecting this sector to the perturba- 
tive high energy regime requires non-perturbative control over the theory, which can 
be achieved by lattice gauge theories. In these lectures the emphasis is on a non- 
perturbative treatment of the theory including its renormalization. Connecting the 
perturbative and the non-perturbative regions is one of the main subjects (II.)- 

Heavy quarks require special care on a lattice with spacing a, simply because their 
mass is of the order of the cutoff, a -1 , or higher. Effective theories can be used in 
this situation, in particular Heavy Quark Effective Theory (HQET) is appropriate for 
hadrons with a single heavy quark. It allows to compute the expansion of their properties 
in terms of the inverse quark mass (in practice the b-quark mass) . The renormalization 
of this effective theory introduces couplings whose number grows with the order of the 
expansion. In order to preserve the predictivity of the theory, these couplings ought to 
be determined from the underlying theory, QCD. Again this step can be seen as the 
renormalization of the effective theory. As will be explained, non-perturbative precision 
is required if one wants to be able to take the continuum limit of the lattice effective 
field theory. Actually it is a general property of the expansion that a 1 jm^ correction 
is only defined once all parts including the matching are done non-perturbatively. 

Renormalization is an ultraviolet phenomenon with relevant momentum scales of 
order a" 1 . Since the QCD coupling becomes weak in the ultraviolet, one may expect to 
be able to perform renormalizations perturbatively, i.e. computed in a power series in 
the bare coupling g§ as one approaches the continuum limit a — ► 0. 2 However, one has 
to care about the following point. In order to keep the numerical effort of a simulation 
tractable, the number of degrees of freedom in the simulation may not be excessively 
large. This means that the lattice spacing a can not be taken very much smaller than the 
relevant physical length scales of the observable that is considered. Consequently the 
momentum scale a -1 that is relevant for the renormalization is not always large enough 
to justify the truncation of the perturbative series. In particular one has to remember 
that the bare coupling vanishes only logarithmically as a — ► 0: ~ l/log(aAQCD)- 
In order to obtain a truly non-perturbative answer, the renormalizations have to be 
performed non-perturbatively. 

Depending on the observable, the necessary renormalizations are of different nature. 
I will use this introduction to point out the different types, and in particular explain the 
problem that occurs in a non-perturbative treatment of scale dependent renormalization. 

2 For simplicity we ignore here the cases of mixing of a given operator with operators of lower dimen- 
sion where this statement does not hold. 
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1.1 Basic renormalization: hadron spectrum 

The calculation of the hadron spectrum starts by choosing certain values for the bare 
coupling, go, and the bare masses of the quarks in units of the lattice spacing, amo,i- 
The flavor index i assumes values i = u, d, s, c, b for the up, down, charm and bottom 
quarks that are sufficient to describe hadrons of up to a few GeV masses. We ignore the 
problem of simulating the b-quark for the moment, neglect isospin breaking and take 
the light quarks to be degenerate, mo jU = m .d = Wo,l- 

Next, from MC simulations of suitable correlation functions, one computes masses 
of five different hadrons H, e.g. H = p, tt, K, D, B for the proton, the pion and the K-,D- 
and B-mesons, 

am H = am H (go,amo :h am 0)S ,amo :C ,am 0:h ) . (1.1) 

The theory is renormalized by first setting m p = m^ p , where is the experimental 
value of the proton mass. This determines the lattice spacing via 

a = (am p )/m c p xp . (1.2) 

Next one must choose the parameters amo,i such that (1.1) is indeed satisfied with the 
experimental values of the meson masses. Equivalently, one may say that at a given 
value of go one fixes the bare quark masses from the condition 

(am H )/(am p ) = m c * p /m c p xp , # = 7r,K,D,B . (1.3) 

and the bare coupling go then determines the value of the lattice spacing through 
eq. (1.2). 

After this renormalization, namely the elimination of the bare parameters in favor 
of physical observables, the theory is completely defined and predictions can be made. 
E.g. the leptonic decay constant, F n , of the pion can be determined, 

F-k = a~ 1 [aF 7r ][l + 0(a)] . (1.4) 

For the rest of this section, I assume that the bare parameters have been eliminated 
and consider the additional renormalizations of more complicated observables. 

Note. Renormalization as described here is done without any reference to pertur- 
bation theory. One could in principle use the perturbative formula for (gAqcd)(<?o) 
for the renormalization of the bare coupling, where Aqcd denotes the A-parameter of 
the theory (in some scheme). Proceeding in this way, one obtains a further prediction 
namely m p /AQCD but at the price of introducing 0(<7q) errors in the prediction of the 
observables. As mentioned before, such errors decrease very slowly as one performs the 
continuum limit. A better method to compute the A-parameter will be discussed later. 
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1.2 Scale dependent renormalization and fundamental parameters of QCD 



As we take the relevant length scales in correlation functions to be small or take the 
energy scale in scattering processes to be high, QCD is better and better approximated 
by weakly coupled quarks and gluons. The strength of the interaction may be measured 
for instance by the ratio of the production rate of three jets to the rate for two jets in 
high energy e + e~ collisions 3 



We observe the following points. 

• The perturbative renormalization group tells us that decreases logarithmi- 
cally with growing energy [i. In other words the renormalization from the bare 
coupling to a renormalized one is logarithmically scale dependent. 

• Different definitions of a are possible; but with increasing energy, a depends less 
and less on the definition (or the process). 

• In a similar way, one may define running quark masses m from combinations of 
observables at high energies. 

• Using a suitable definition (scheme), the //-dependence of a and m can be de- 
termined non-perturbatively and at high energies the short distance parameters 
a and fn can be converted to the renormalization group invariants using pertur- 
bation theory in a. Being defined non-perturbatively, the latter are the natural 
fundamental parameters of QCD. 

Explaining these points in detail is the main objective of the second lecture. 
1.3 Irrelevant operators 

Another category of renormalization is associated with the removal of lattice discretiza- 
tion errors such as the linear a-term in eq. (1.4). Following Symanzik's improvement 
program, this can be achieved order by order in the lattice spacing by adding irrelevant 
operators, i.e. operators of dimension larger than four, to the lattice Lagrangian [1]. 
The coefficients of these operators are easily determined at tree level of perturbation 
theory, but in general they need to be renormalized. We will explain the general idea of 
the non-perturbative determination of the coefficients arising at order a and then briefly 
review the present status of 0(a) improvement. 

Note also the alternative approach of removing lattice artifacts order by order in 
the coupling constant but non-perturbatively in the lattice spacing a described in the 

3 One should really use some rather inclusive process, e.g. one computable directly in the Euclidean 
theory. For explaining the principle we ignore this issue. 
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lectures by Peter Hasenfratz. Linear effects in a are automatically absent if the lattice 
regularization has enough chiral symmetry. Indeed chiral symmetry can be kept exactly 
in the discretized theory [2-4], but these theories are rather expensive to simulate. On 
the other hand also the "twisted mass" regularization [5, 6] is automatically 4 0(a)- 
improved [7] (see the appendix of [8] and Stefan Sint's lectures at this school for a 
simple argument), but at the price of the violation of isospin symmetry. 

1.4 Heavy Quark Effective Theory 

This theory is very promising for B-physics. It approximates heavy-light bound state 
properties systematically in an expansion of AQco/^b) a small expansion parameter. 
A non-trivial issue is the renormalization of the theory. Already at the lowest order in 
1/mb, the associated uncertainties are significant if renormalization is treated pertur- 
batively. At that order renormalization can be carried out by the methods discussed in 
the second lecture [9-11], but when one includes 0(l/mb) corrections one has to deal 
in addition with the mixing of operators of different dimensions. 5 The continuum limit 
of the effective theory then exists only if the power divergent mixing coefficients are 
computed non-perturbatively. 

In the third lecture we will explain these issues in detail. We will formulate HQET 
non-perturbatively. The power divergent mixing coefficients can then be determined by 
matching the theory to QCD. A possible strategy will be explained. As an example we 
will show the computation of the b-quark mass including corrections. 



4 "Automatically" still means that the standard mass term has to be tuned to zero, but that can be 
done by the use of the PCAC relation. 

Note that the computation of an order Aqcd term in the renormalized quark mass already consti- 
tutes a 0(l/m,b) correction to the leading term, although it is done in static approximation. 
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Lecture I. 

The Schrodinger functional and O(a)-improvement of lattice QCD 



1.1 The Schrodinger functional (SF) 

For various applications, for instance scale dependent renormalization in QCD, 0(a)- 
improvement and Heavy Quark Effective Theory, we need QCD in a finite volume with 
boundary conditions suitable for (easy) perturbative calculations and MC simulations. 
These are provided by the SF of QCD, which we introduce below. For a while we 
restrict the discussion to the pure gauge theory. In this part the presentation follows 
closely [12]; we refer to this work for further details as well as proofs of the properties 
described below. 
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space 

(LxLxL box with periodic b.c.) 

Figure 1: Illustration of the Schrodinger functional. 



1.1.1 Definition 

Here, we give a formal definition of the SF in the Yang-Mills theory in continuum space- 
time, noting that a rigorous treatment is possible in the lattice regularized theory. 

Space-time is taken to be a cylinder illustrated in Fig. 1. We impose Dirichlet 
boundary conditions for the vector potentials 6 in time, 

= {is :=° • ^ 

where C, C are classical gauge potentials and A A denotes the gauge transform of A, 
A£(x) = A(x)^(x)A(x)- 1 + A(x)a fc A(x)- 1 , A £ SU(jV) . (1.1.2) 



J We use anti-hermitian vector potentials. 
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In space, we impose periodic boundary conditions, 

A k (x + Lk) = A k (x), A(x + Lk) = A(x) . (1.1.3) 
The (Euclidean) partition function with these boundary conditions defines the SF, 

Z[C',C] = Jb[A] jB[A]e~ SG[A] , (1.1.4) 

S G [A] = f d 4 xtr {F^F^}, 

Ffj,u = d^A u — dyA^ + [Ap, A u \ , 
B[A] = J] <UJ(s), D[A]=ndA(x) . 

x,fi,a x 

Here dA(x) denotes the Haar measure of SU(iV). It is easy to show that the SF is a 
gauge invariant functional of the boundary fields, 

Z[C' n ',C n ] = Z[C',C] , (1.1.5) 

where also large gauge transformations are permitted. The invariance under the latter is 
an automatic property of the SF defined on a lattice, while in the continuum formulation 
it is enforced by the integral over A in eq. (1.1.4). 

1.1.2 Quantum mechanical interpretation 

The SF is the quantum mechanical transition amplitude from a state \C) to a state \C) 

after a (Euclidean) time L. To explain the meaning of this statement, we introduce 

the Schrodinger representation. The Hilbert space consists of wave-functionals ^f[A] 

which are functionals of the spatial components of the vector potentials, ^4^(x). The 

canonically conjugate field variables are represented by functional derivatives, ^(x) = 
1 6 



and a scalar product is given by 



= J D[A]*[A]**'[A] ) B[A] = l[ dA£(x) . (1.1.6) 



x,fc,a 



The Hamilton operator, 



H = J o d 3 x||^(x)^(x) + i ^F^(x)^(x)| , (1.1.7) 

commutes with the projector, IP, onto the physical subspace of the Hilbert space (i.e. 
the space of gauge invariant states), where IP acts as 

TPif>[A] = f B[A}^j[A a ] . (1.1.8) 
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Finally, each classical gauge field defines a state \C) through 

(C|¥) = *[C] . (1.1.9) 
After these definitions, the quantum mechanical representation of the SF is given by 

Z[C',C] = (C'\e~ HT W>\C) 

oo 

= J2 e ~ EnTq, n[C']* n [C]* . (1.1.10) 
n=0 

In Wilson's original lattice formulation, eq. (1. 1.10) can be derived rigorously and is 
valid with real energy eigenvalues E n . 

1.1.3 Background field 

A complementary aspect of the SF is that it allows a treatment of QCD in a color 
background field in an unambiguous way. Let us assume that we have a solution B of 
the equations of motion, which satisfies also the boundary conditions eq. (1. 1.1). If, in 
addition, 

S[A] > S[B] (1.1.11) 

for all gauge fields A that are not equal to a gauge transform B of B, then we call 
B the background field (induced by the boundary conditions). Here, O(x) is a gauge 
transformation defined for all x in the cylinder and its boundary and B^ 1 is the corre- 
sponding generalization of eq. (1.1.2). Background fields B, satisfying these conditions 
are known; we will describe a particular family of fields, later. 

Due to eq. (1. 1.11), fields close to B dominate the path integral for weak coupling 
go and the effective action, 

T[B] = -]nZ[C',C] , (1.1.12) 
has a regular perturbative expansion, 

T[B] = \r [B}+T 1 [B}+g 2 T 2 [B} + ... , (1.1.13) 
% 

T [B] = g 2 S[B] . 

Above we have used that due to our assumptions, the background field, B, and the 
boundary values C, C are in one-to-one correspondence and have taken B as the argu- 
ment of r. 
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1.1.4 Perturbative expansion 

For the construction of the SF as a renormalization scheme, one needs to study the 
renormalization properties of the functional Z. Liischer, Narayanan, Weisz and Wolff 
have performed a one-loop calculation for arbitrary background field [12]. The calcula- 
tion is done in dimensional regularization with 3 — 2e space dimensions and one time 
dimension. One expands the field A in terms of the background field and a fluctuation 
field, q, as 

A^x) = B^x) + g q^x) . (1.1.14) 

Then one adds a gauge fixing term ("background field gauge") and the corresponding 
Fadeev-Popov term. Of course, care must be taken about the proper boundary condi- 
tions in all these expressions. Integration over the quantum field and the ghost fields 
then gives 

ri[J3] = ilndetAi -lndetA , (1.1.15) 

where Ax is the fluctuation operator and Ao the Fadeev-Popov operator defined in [12]. 
The result can be cast in the form 

Ti[B] = -^r [B]+O(l) , (1.1.16) 

with the important result that the only (for s — ► 0) singular term is proportional to To- 
After renormalization of the coupling, i.e. the replacement of the bare coupling by 

9ms via 

9l= P 2£ 9mW + ^)s^s^)], zi(e)=-j, (1.1.17) 
the effective action is finite, 

T[B] e=0 = jj--&o 

-K'(0|Ai) 
C'(0|A) = A C(s | A) 

dfi S : 

Here, £'(0|A) is a complicated functional of B, which is not known analytically but can 
be evaluated numerically for specific choices of B. 

The important result of this calculation is that (apart from field independent terms 
that have been dropped everywhere) the SF is finite after eliminating qq in favor of 
</jyjg. The presence of the boundaries does not introduce any extra divergences. In the 
following section we argue that this property is correct in general, not just in one-loop 
approximation . 



hV 



l 

16^ 



To [B] 



+ C'(0|A ) + O(5|s) 
C(s|A) = TrA~ 



(1.1.18) 



8 



1.1.5 General renormalization properties 

The relevant question here is whether local quantum field theories formulated on space- 
time manifolds with boundaries develop divergences that are not present in the absence 
of boundaries (periodic boundary conditions or infinite space-time). In general the 
answer is "yes, such additional divergences exist". In particular, Symanzik studied 
the </> 4 -theory with Schrodinger functional boundary conditions [13,14]. He presented 
arguments that to all orders of perturbation theory the Schrodinger functional is finite 
after 

• renormalization of the self-coupling, A, and the mass, m, 

• and the addition of the boundary counter-terms 



In addition to the standard renormalizations, one has to add counter-terms formed 
by local composite fields integrated over the boundaries. One expects that in general, 
all fields with dimension d < 3 have to be taken into account. Already Symanzik 
conjectured that counter-terms with this property are sufficient to renormalize the SF 
of any quantum field theory which is renormalizable when no boundaries are present. 

Since this conjecture forms the basis for many applications of the SF, we note a 
few points concerning its status. 

• As mentioned, a proof to all orders of perturbation theory does not exist. An 
application of power counting in momentum space in order to prove the conjecture 
is not possible due to the missing translation invariance. 

• There is no gauge invariant local field with d < 3 in the Yang-Mills theory. 
Consequently no additional counter-term is necessary in accordance with the 1- 
loop result described in the previous subsection. 

• In QCD it has been checked also by explicit 2-loop calculations [15, 16]. MC 
simulations give further support for its validity beyond perturbation theory; we 
give examples in the second lecture. 

Although a general proof is missing, there is little doubt that Symanzik's conjecture 
is valid in general. Concerning QCD, this puts us into the position to give an elegant 
definition of a renormalized coupling in finite volume. 

1.1.6 Renormalized coupling 

For the definition of a running coupling we need a quantity which depends only on 
one scale. We choose LB such that it depends only on one dimensionless variable rj. 
In other words, the strength of the field is scaled as 1/L. The background field is 




(1.1.19) 
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assumed to fulfill the requirements of Sect. 1.1.3. Then, following the above discussion, 
the derivative 



r'[B] = —T[B] , 



(1.1.20) 



is finite when it is expressed in terms of a renormalized coupling like g^jg but V is 
defined non-perturbatively. From eq. (1.1.13) we read off immediately that a properly 
normalized coupling is given by 



g 2 (L) = T' [B}/T'[B] . 



(1.1.21) 



Since there is only one length scale L, it is evident that g defined in this way runs with 
L. 

A specific choice for the gauge group SU(3) is the abelian background field induced 
by the boundary values [17] 
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(1.1.22) 
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(1.1.23) 



In this case, the derivatives with respect to rj are to be evaluated at i] = 0. The 
associated background field, 

B = 0, B k =[x C' k + (L-x )C k ]/L, fc = l,2,3, (1.1.24) 

has a field tensor with non- vanishing components 

G ok = d B k = (C' k - C k )/L, = 1,2,3 . (1.1.25) 

It is a constant color-electric field. 



1.1.7 Quarks 

In the end, the real interest is in the renormalization of QCD and we need to consider 
the SF with quarks. We restrict our discussion to the original formulation of S. Sint [18]. 

Special care has to be taken in formulating the Dirichlet boundary conditions for 
the quark fields; since the Dirac operator is a first order differential operator, the Dirac 
equation has a unique solution when one half of the components of the fermion fields are 
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specified on the boundaries. Indeed, a detailed investigation shows that the boundary 
condition 

P+^L =0 = p, P-ip\ X0 =L = p' , P± = i(l± 7 o), (1-1.26) 
^-P-U- =o = P, 4>P+\x =L = P, (1.1.27) 

lead to a quantum mechanical interpretation analogous to eq. (1. 1.10). The SF 

Z[C',p',p';C,p,p] = J D[^]D[^]D[^]e- 5 [ A ^l (1.1.28) 

involves an integration over all fields with the specified boundary values. The full action 
may be written as 

S[A^^\ = S G ®,1>] + S F [Atf,1>] 

S F = /Vx^(x)[ 7/x I^ + m]#c) (1.1.29) 

- J d 3 x $(x)P^(x)] xo = - J d 3 x $(x)P + i;(x)} X0=L , 

with Sq as given in eq. (1.1.4). In eq. (1.1.29) we use standard Euclidean 7-matrices. 
The covariant derivative, D^, acts as D^ip(x) = d^ip^x) + A fl (x)ijj(x). 

Let us now discuss the renormalization of the SF with quarks. In contrast to the 
pure Yang-Mills theory, gauge invariant composite fields of dimension three are present 
in QCD. Taking into account the boundary conditions one finds [18] that the counter- 
terms, 

^P_VUo=o and ^P + i>\ X0=L , (1.1.30) 

have to be added to the action with weight 1 — to obtain a finite renormalized 
functional. These counter-terms are equivalent to a multiplicative renormalization of 
the boundary values, 

PR = 1/2 Pi ■■■ » PR = Z b 1/2 P ■ (1.1.31) 

It follows that - apart from the renormalization of the coupling and the quark mass 
- no additional renormalization of the SF is necessary for vanishing boundary values 
p, . . . , p' . So, after imposing homogeneous boundary conditions for the fermion fields, a 
renormalized coupling may be defined as in the previous subsection. 

As an important aside, we point out that the boundary conditions for the fermions 
introduce a gap into the spectrum of the Dirac operator (at least for weak couplings). 
One may hence simulate the lattice SF for vanishing physical quark masses. It is then 
convenient to supplement the definition of the renormalized coupling by the requirement 
m = 0. In this way, one defines a mass-independent renormalization scheme with simple 
renormalization group equations. In particular, the /3-function remains independent of 
the quark mass. 
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1.1.7.1 Correlation functions 

are given in terms of the expectation values of any product O of fields, 

(O) = {\\ B[A]B[i;}B[^]Oe- s ^^} , (1.1.32) 

I £ J ) p'=p'=p=p=0 

evaluated for vanishing boundary values p, . . . ,p'. Apart from the gauge field and the 
quark and anti-quark fields integrated over, O may involve the "boundary fields" [19] 

5 * 
C(x) = — -, C(x 



«*> = ski' ?(x) = ski ■ (LL33) 

An application of fermionic correlation functions including the boundary fields is the 
definition of the renormalized quark mass in the SF scheme to be discussed next. 

1.1.7.2 Renormalized mass 

Just as in the case of the coupling constant, there is a great freedom in defining renor- 
malized quark masses. A natural starting point is the PCAC relation which expresses 
the divergence of the axial current 7 , 

Al(x) = ?(x) 7 57„KV(x) , (1-1.34) 

(for simplicity we have chosen just Nf = 2 degenerate flavors and the Pauli matrix r a 
acts in this flavor space), in terms of the associated pseudo-scalar density, 

P a (x) = ?(x) 7 4t>(x) , (1.1.35) 

via 

dpA^x) = 2mP a (x) . (1.1.36) 

This operator identity is easily derived at the classical level (cf. Sect. 1.2). After 
renormalizing the operators, 

(A R )« = Z A Al , 

= Z P P a , (1.1.37) 

a renormalized current quark mass may be defined by 

TiR = fn = -J~fn . (1.1.38) 
Z P 



7 The reader is not to confuse A^(x), with the gauge vector potentials A^(x) 
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Here, m, is to be taken from eq. (1.1.36) inserted into an arbitrary correlation function 
and Za can be determined from a proper chiral Ward identity [20-22]. Note that m 
does not depend on which correlation function is used because the PCAC relation is an 
operator identity. The definition of fn is completed by supplementing eq. (1.1.37) with 
a specific normalization condition for the pseudo-scalar density. The running mass fn 
then inherits its scheme- and scale-dependence (n) from the corresponding dependence of 
Pr. Such a normalization condition may be imposed through infinite volume correlation 
functions. Since we want to be able to compute m(/i) for large energy scales fi, we do, 
however, need a finite volume definition (see Sect. II. 1.2). This is readily given in terms 
of correlation functions in the SF. 




space 



Figure 2: /p (left) and f\ (right) in terms of quark propagators. 
To start with, let us define (iso-vector) pseudo-scalar fields at the boundary of the 

SF, 

O a = J d 3 u J d 3 vC(u) 75 ^r a C(v), 

O la = |d 3 uy>vC'(u)T5ir°C'(v) , (L1.39) 

to build up the correlation functions 

/p(x ) = -\{P a {x)O a ) , 

f x = (0 ,a O a ) , (1.1.40) 

which are illustrated in Fig. 2. 
We then form the ratio 

Z P = const. v7i//p(a:)L=L/2 » (1.1.41) 

such that the renormalization of the boundary quark fields, eq. (1.1.31), cancels out. 
The proportionality constant is to be chosen such that Zp = 1 at tree level. To define 
the scheme completely one needs to further specify the boundary values C, C and the 
boundary conditions for the quark fields in space. These details are of no importance, 
here. 
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We rather mention some more basic points about this renormalization scheme. 
Just like in the case of the running coupling, the only physical scale that exists in our 
definitions eq. (1.1.38), eq. (1.1.41) is the linear dimension of the SF, the length scale, 
L. So the mass rn(fi) runs with \x = 1/L. We have already emphasized that g is to be 
evaluated at zero quark mass. It is advantageous to do the same for Z-p. In this way we 
define a mass-independent renormalization scheme, with simple renormalization group 
equations. 

By construction, the SF scheme is non-perturbative and independent of a specific 
regularization. For a concrete non-perturbative computation, we do, however, need to 
evaluate the expectation values by a MC-simulation of the corresponding lattice theory. 
We proceed to introduce the lattice formulation of the SF. 

1.1.8 Lattice formulation 

A detailed knowledge of the form of the lattice action is not required for an under- 
standing of the following sections. Nevertheless, we give a definition of the SF in lattice 
regularization. This is done both for completeness and because it allows us to obtain a 
first impression about the size of discretization errors. 

We choose a hyper-cubic Euclidean lattice with spacing a. A gauge field U on 
the lattice is an assignment of a matrix U(x,fi) £ SU(iV) to every lattice point x and 
direction fx = 0, 1, 2, 3. Quark and anti-quark fields, ip{x) and ip(x), reside on the lattice 
sites and carry Dirac, color and flavor indices as in the continuum. To be able to write 
the quark action in an elegant form it is useful to extend the fields, initially defined only 
inside the SF manifold (cf. Fig. 1) to all times xo by "padding" with zeros. In the case 
of the quark field one sets 

ijj{x) = if xo < or xo > L, 

and 

P-^(x)\x o =0 = P+1p(x)\x =L = 0, 

and similarly for the anti-quark field. Gauge field variables that reside outside the 
manifold are set to 1. 

We may then write the fermionic action as a sum over all space-time points without 
restrictions for the time-coordinate, 

S ¥ [U^,i>] = a 4 ^(£> + m )V, (1-1-42) 

X 

and with the standard Wilson-Dirac operator, 

D = \ E{7m(V* + V„) - aV* V M } . (1.1.43) 

/j,=0 
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Here, forward and backward covariant derivatives, 



V^(z) = \[U{x,^{x + a(i)-^{x)\, (1.1.44) 
V* ip(x) = \[^{x) - U(x - aft, fi)- 1 ^ - aft)] , (1.1.45) 

are used and tuq is to be understood as a diagonal matrix in flavor space with elements 
m J . 

The gauge field action Sq is a sum over all oriented plaquettes p on the lattice, 
with the weight factors w(p), and the parallel transporters U(p) around p, 

S g[U] = W 1 " U(P)} ■ (1.1.46) 

% P 

The weights w{p) are 1 for plaquettes in the interior and 

w (p\ — 1 if P is a spatial plaquette at xq = or xq = L, _ ^ 

I ct if p is timedike and attached to a boundary plane. 

The choice c s = ct = 1 corresponds to the standard Wilson action. However, these 
parameters can be tuned in order to reduce lattice artifacts, as will be briefly discussed 
below. 

With these ingredients, the path integral representation of the Schrodinger func- 
tional reads [18], 

Z = jD[ifj]D[4>]D[U]e- s , S = S F + S G , (1.1.48) 
B[U]=HdU(x,n) , 

X,fl 

with the Haar measure dll. 



1.1.8.1 Boundary conditions and the background field. 

The boundary conditions for the lattice gauge fields may be obtained from the continuum 
boundary values by forming the appropriate parallel transporters from x + ak to x at 
xq = and xq = L. For the constant abelian boundary fields C and C that we 
considered before, they are simply 

U(x, fc) | xo =o = exp(aC fc ), U(x, k)\ X0=L = exp(aC£), (1.1.49) 

for k = 1, 2, 3. All other boundary conditions are as in the continuum. 

For the case of eq. (I.1.22),eq. (1.1.23), the boundary conditions (1.1.49) lead to 
a unique (up to gauge transformations) minimal action configuration V, the lattice 
background field. It can be expressed in terms of B (1.1.24), 

V(x,fi) = exp{aB fi (x)} . (1.1.50) 
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1.1.8.2 Lattice artifacts. 

Now we want to get a first impression about the dependence of the lattice SF on the 
value of the lattice spacing. In other words we study lattice artifacts. At lowest order 
in the bare coupling we have, just like in the continuum, 



1 



ro[^ + 0((<7o) ), T [V]^g 2 S G [V] 



(1.1.51) 



Furthermore one easily finds the action for small lattice spacings, 



Sg[V] 



1 + 

o N 

% Q = l 



3L 



4 N 



sm 



2L 2 



l + (l-c t )^ + 0(a 4 ) 



(1.1.52) 



We observe: at tree-level of perturbation theory, all linear lattice artifacts are removed 
when one sets ct = 1. Beyond tree-level (and in the theory without quarks), one has to 
tune the coefficient Ct as a function of the bare coupling. We will show the effect, when 
this is done to first order in g$, below. Note that the existence of linear 0(a) errors in 
the Yang- Mills theory is special to the SF; they originate from dimension four operators 
^bfc-^bfc an d FkiFki which are irrelevant terms (i.e. they carry an explicit factor of the 
lattice spacing) when they are integrated over the surfaces. c s , which can be tuned 
to cancel the effects of F^iF^i, does not appear for the electric field that we discussed 
above. 

Once quark fields are present, there are more irrelevant operators that can gener- 
ate 0(a) effects as discussed in detail in [19]. Here we emphasize a different feature of 
eq. (1.1.52): once the 0(a)-terms are canceled, the remaining a-effects are tiny. This 
special feature of the abelian background field is most welcome for the numerical com- 
putation of the running coupling; it allows for reliable extrapolations to the continuum 
limit. 



1.1.8.3 Explicit expression for V . 

Let us finally explain that V is an observable that can easily be calculated in a MC 
simulation. From its definition we find immediately 



r' 



/ D[ip]D[il>]D[U]e- 



/dS G \ /dS F \ 



(1.1.53) 



drj [J l7-j ly- j j - j \ Qrj / \ drj / 

The derivative evaluates to the (color 8 component of the) electric field at the 
boundary, 

dS a 2 



■a^{^(x)-^)'( X )} , 



dl] g$L 

£|(x) = -^Re tr {i\ fi U{x,k)U{x + ak,Q)U{x + aQ 1 k)- l U{x,{))- 1 \ 



(1.1.54) 
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where As = diag(l, —1/2, —1/2). (A similar expression holds for (_E|)'(x)). The second 
term ^jp, which is only non-zero in the 0(a)-improved formulation is numerically less 
relevant. An explicit expression is given in [23]. 

The renormalized coupling is related to the expectation value of a local operator; 
no correlation function is involved. This means that it is easy and fast in computer time 
to evaluate it. It further turns out that a good statistical precision is reached with a 
moderate size statistical ensemble. 

1.1.9 More literature 

We here give some guide for further reading on the SF. Independently of the work of 
Symanzik, G. C. Rossi and M. Testa discussed different boundary conditions imposed 
at fixed time [24,25]. The renormalization properties of that functional have not yet 
been discussed. 

There are also rather recent developments. Different formulations of the lattice 
Schrodinger functional with overlap fermions satisfying the Ginsparg Wilson relation 
have been found by Y. Taniguchi using an orbifold construction [26], and by M. Liischer 
using a general universality argument concerning QFT's with boundaries [14]; see also 
[27]. As the Schrodinger functional breaks chiral symmetry by the boundary conditions, 
it is relevant into which direction in flavor space the mass term is introduced. The 
Schrodinger functional with a twisted mass term [5, 6] and the boundary conditions 
specified above differs from the SF with a standard mass term (at finite quark mass). 
S. Sint found a modification of the boundary conditions, which yields the standard SF 
as the continuum limit of the lattice theory with a twisted mass [28]. An even number 
of flavors is required in this formulation. It offers also advantages in the massless limit, 
where "automatic bulk 0(a)-improvement" is achieved after the tuning of one counter- 
term. It is discussed in detail in S. Sint's lectures at this school. Another Schrodinger 
functional with automatic bulk 0(a)-improvement is proposed in [29]. 

1.2 Chiral symmetry and O(a)-improvement 

The main focus of this section is on the 0(a) improvement of Wilson's lattice QCD. 
However, we also mention the finite normalization of isovector currents. Both of these 
problems have the same origin, namely that chiral symmetry is broken in Wilson's 
regularization and then also the same solution: chiral Ward identities. The possibility 
to use these to normalize the currents has first been discussed by Refs. [20,30]. Here, we 
describe their application in the computation of the 0(a)-improved action and currents. 
A difference to the aforementioned work is that we insist that only on-shell improvable 
correlation functions are used in the normalization conditions in order to be compatible 
with on-shell improvement. 

Before going into more details, we would like to convey the general idea of the 
application of chiral Ward identities. For simplicity we assume an isospin doublet of 
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mass-degenerate quarks. Consider first a regularization of QCD which preserves the 
full SU(2)v xSU(2)a flavor symmetry as it is present in the continuum Lagrangian of 
massless QCD. Indeed, such regularizations exist [2-4], see Peter Hasenfratz' lectures. 

We can derive chiral Ward identities in the Euclidean formulation of the theory. 
These then provide exact relations between different correlation functions. Immediate 
consequences of these relations are that there are currents V^, which do not get 
renormalized (Z\ = Zy = 1) and the quark mass does not have an additive renormal- 
ization. 

In a general discretization, such as the Wilson formulation, lattice QCD does not 
have the SU(2)a flavor symmetry for finite values of the lattice spacing. Then, the Ward 
identities are not satisfied exactly. From universality, we do, however, expect that the 
correlation functions may be renormalized such that they obey the same Ward identities 
as before - up to 0(a) corrections. Therefore we may impose those Ward identities for 
the renormalized currents, to fix their normalizations. 

Furthermore, following Symanzik [1], it suffices to a add a few local irrelevant terms 
to the action and to the currents in order to obtain an improved lattice theory, where 
the continuum limit is approached with corrections of order a 2 . The coefficients of these 
terms can be determined by imposing improvement conditions. For example one may 
require certain chiral Ward identities to be valid at finite lattice spacing a. 

1.2.1 Chiral Ward identities 

For the moment we do not pay attention to a regularization of the theory and derive 
the Ward identities in a formal way. As mentioned above these identities are exact in a 
regularization that preserves chiral symmetry. To derive the Ward identities, one starts 
from the path integral representation of a correlation function and performs the change 
of integration variables 

^( x ) _> e i2 ? [e A(^+4M]^( x ) 

= ^(x)+ie a A (x)5 a A ^(x)+iey(x)5yiP(x) , 

= ^(x) + ie a A (x)5 a A ^(x) + iey(x)Sy^(x) , (1.2.1) 
where we have taken e A (x),€y(x) infinitesimal and introduced the variations 

5 a A ^{x) = \T a lb ^(x), 5 a ^(x) =^(x) 75 ±T a . (1.2.2) 

The Ward identities then follow from the invariance of the path integral representation 
of correlation functions with respect to such changes of integration variables. They 
obtain contributions from the variation of the action and the variations of the fields in 
the correlation functions. The variations of the currents 

V;{x) =^{x) lfl \r a ^(x) (1.2.3) 
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and A^, eq. (1.1.34), is given by 



-ie abc VC(x), 
-ie^AUx), 



-u 



-ic 



bc V'(x) . 



(1.2.4) 



They form a closed algebra under these variations. 

Since this is convenient for our applications, we 
write the Ward identities in an integrated form. 
Let R be a space-time region with smooth bound- 
ary dR. Suppose Oi nt and O cxt are polynomials in 
the basic fields localized in the interior and exte- 
rior of R respectively (see left). The general vector 
current Ward identity then reads 




ext 



d^(ar) (V?(x)O int O, 



dR 



while for the axial current one obtains 



ext 



((^o int ) a 



ext/ 



(1.2.5) 



J ^(x) (A°(x)O int O cxt ) = - {(5%O int ) O e xt) (1-2.6) 

+2m / d 4 x (P a (x)O int O cxt ) . 
Jr 

Here volume integrals over for example d^A^x) have been changed to surface integrals. 
The integration measure da^(x) points along the outward normal to the surface dR and 
P a (x) was defined in eq. (1.1.35). 

We may also write down the precise meaning of the PCAC-relation eq. (1.1.36). It 
is eq. (1.2.6) in a differential form, 



d l ,A a {x)-2mP a {x) O, 



'ext 







(1.2.7) 



where now ex t may have support everywhere except for at the point x. 

Going through the same derivation in the lattice regularization, one finds equations 
of essentially the same form as the ones given above, but with additional terms [20]. At 
the classical level these terms are of order a. More precisely, in eq. (1-2.7) the important 
additional term originates from the variation of the Wilson term, a 5\{ip\/*yV ^ip) , and 
is a local field of dimension 5. Such 0(a)-corrections are present in any observable 
computed on the lattice and are no reason for concern. However, as is well known 
in field theory, such operators mix with the ones of lower and equal dimensions when 
one goes beyond the classical approximation. In the present case, the dimension five 
operator mixes among others also with d^A^(x) and P a (x). This means that part of the 
classical 0(a)-terms turn into 0(<7q) hi the quantum theory. The essential observation 
is now that this mixing can simply be written in the form of a renormalization of the 
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terms that are already present in the Ward identities, since all dimension three and four 
operators with the right quantum number are already there. 

We conclude that the identities, which we derived above in a formal manner, are 
valid in any proper lattice regularization after 

• replacing the bare fields A, V, P and quark mass mo by renormalized ones, where 
one must allow for the most general renormalizations, 

(A R )« = Z A A; , (Va)« = ZyVf , 

(PrT = Z P P a , tor = Z m m q , TO q = to - m c , 

• allowing for the usual 0(a) lattice artifacts. 

Note that the additive quark mass renormalization m c diverges like 0(g 2 /a) for dimen- 
sional reasons. 

As a result of this discussion, the formal Ward identities may be used to determine 
the normalizations of the currents. We refer the reader to [21,31] for details and explain 
here the general idea how one can use the Ward identities to determine improvement 
coefficients. 

1.2.2 On-shell O(a)-improvement 

1.2.2.1 Motivation 

Let us first recall why one wants to remove lattice spacing effects linear in a. The prime 
reason is as follows. If linear effects are present, one has to vary a in the numerical 
simulations over a large range in order to be able to get a reasonable estimate of their 
magnitude. 8 In contrast if the cutoff effects are 0(a 2 ), a range of 0.05 fm < a < 0.1 fm 
typically allows to check well whether they contribute significantly. In fact a reasonably 
well controlled extrapolation to the continuum can then be done allowing for a term 
proportional to a 2 and also a smaller range in a may be sufficient. Examples can be 
found e.g. in [33]. In addition, it does turn out that linear a effects can be quite large. 
Let us give here just two examples. 

The first one is the current quark mass to defined by the PCAC relation. As 
detailed below, its value is independent of kinematical variables such as the boundary 
conditions. Dependences on such variables are pure lattice artifacts. We examined the 
current quark mass in the valence approximation by numerical Monte Carlo simulations 
and found large lattice artifacts even for quite small lattice spacings [32] (cf. fig. 3). 

The second example is the mass of the vector meson, made dimensionless by multi- 
plying with ro . This quantity has large cutoff effects in the quenched approximation [34] . 
Depending on the quark mass, a-effects of around 20% and more are seen at a ~ 0.15 fm; 
see for example Fig. 1 of [34], Fig. 1 of [35]. 

8 Obviously it does not really help to have a large range by considering large values of a. Then 
one enters the regime where either the higher order terms are significant or - more likely - the whole 
asymptotic expansion in a breaks down. 
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Figure 3: Dependence of current quark mass m on the boundary condition and the time 
coordinate [32]. The calculation is done in the quenched approximation on a 16 x 8 3 
lattice at (3 = 6.4, which corresponds to a lattice spacing of a ~ 0.05 fm. "Boundary 
values" refer to the gauge field boundary conditions in the SF. Their values are given 
in [32]. 

1.2.2.2 A warning from two dimensions 

Before entering the discussion of the O(a) improvement programme, we mention some 
unexpected results from thorough examinations of 2-d O(N) sigma models. The the- 
oretical basis for the discussion and removal of a-effects is Symanzik's effective theory, 
see Sect. 1.2.2.3. O(N) sigma models were the second class of models investigated by 
Symanzik in order to establish this theory. For these models the basic statement is that 
(up to logarithmic modifications) the cutoff effects are quadratic in a, when a is small 
enough. 

It therefore came as a surprise that Hasenfratz and Niedermayer found in a nu- 
merical study of the Liischer-Weisz- Wolff (LWW) renormalized coupling [37] of the 2-d 
0(3) sigma model that its step scaling function shows an a-dependence which is roughly 
linear in a for quite small a (large correlation length) [38]. With a further improved 
algorithm, a Bern-Berlin collaboration confirmed this behavior with even higher preci- 
sion and smaller lattice spacing [36] . We show their result in Fig. 4. On the other hand 
it was known that the cutoff effects of the step scaling function are 0(a 2 ) in the large 
N limit of the O(N) models [37]. Subsequent numerical studies for N = 4, 8 showed 
no conclusive results: just like in the N = 3 case, the cutoff effects look linear when 
judged by eye, but they can also be fitted with 0(a 2 ) functions, in particular when the 
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Figure 4: Lattice spacing dependence of the step scaling function of the LWW coupling 
in the 2-d 0(3) sigma model [36] for coupling uq = 1.0595. The data points with the 
smaller cutoff effects are for the standard nearest neighbor action. 

expected logarithmic modifications are taken into account [39]. 

Later the 1/N correction was worked out at finite lattice spacing [39]. Recall that 
at order (1/N)° one has 0(a 2 ) effects. The cutoff effect proportional to 1/N is shown in 
Fig. 5. Over a large range in a it is almost a linear function of a, but close to a = it is 
dominated by an 0(a 2 ) term. Thus our personal conclusion is that there is no conflict 
with Symanzik's effective theory in the O(N) models. One should also note that all the 
a-effects discussed here are rather small. 

However, there is a clear warning that, depending on model and observable, the 
lattice spacing may have to be rather small for the leading correction term in the effective 
theory to dominate. On the more practical side, long continuum extrapolations with 
significant slopes may be dangerous, since in QCD we do not have much information 
where the asymptotic expansion in a is accurate [40]. This is one of the reasons why 
we will spend much time on understanding the cutoff effects in the QCD step scaling 
function of the coupling in Sect. II. 2. 2 and Sect. II. 2. 3. 

1.2.2.3 Symanzik's local effective theory (SET) 

In the following explanation of the theory we follow quite closely [19]. We consider QCD 
on an infinitely extended lattice with two degenerate light Wilson quarks of bare mass 
mo [41]. The action is then given as in Sect. 1.1.8 except that no boundary conditions 
or boundary terms are necessary 

Quite some time ago, Symanzik provided arguments that a lattice theory can be 
described in terms of a local effective theory, when the lattice spacing is small enough [1] . 
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Figure 5: Coefficient of 1 /N in the 1/iV-expansion of the cutoff effects of the step scaling 
function of the LWW coupling of 2-d 0(N) sigma models. Graph prepare by U. Wolff 
based on [39]. 

The effective action, 

5 eff = S + aSi + a 2 S 2 + ..., (1.2.8) 

has as a leading order, Sq, the action of the continuum theory 9 . The terms Sk, k = 
1,2,..., are space-time integrals of Lagrangians Ck{x). These are given as general 
linear combinations of local gauge-invariant composite fields which respect the exact 
symmetries of the lattice theory and have canonical dimension 4 + k. We use the 
convention that explicit (non-negative) powers of the quark mass m are included in the 
dimension counting. A possible basis of fields for the Lagrangian C±(x) reads 

Oi = ^cr^F^ip, 

3 = mtr{F^F^}, (1.2.9) 

05 = m 2 ip'^, 

where is the gluon field tensor and a^ u = | [7^,7^] . 

When considering correlation functions of local gauge invariant fields the action is 
not the only source of cutoff effects. If <p(x) denotes such a lattice field (e.g. the axial 
density or the isospin currents), one expects the connected n-point function 

G n (xi, ...,x n ) = (Z^ (<P( X1 ) . . . 0(x n )) con (1.2.10) 

9 If more rigor is desired one may define it on a lattice with spacing e«a. 
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to have a well-defined continuum limit, provided the renormalization constant is 
correctly tuned and the space-time arguments kept at a physical distance 

from each other. 

In the effective theory the renormalized lattice field Z ( p<j)(x) is represented by an 
effective field, 

0cs(x) = <fo{x) + aMx) + a 2 Mx) + ■■■, (1.2.11) 

where the 4>k(x) are linear combinations of composite, local fields with the appropriate 
dimension and symmetries. For example, in the case of the axial current (1.1.34), (pi is 
given as a linear combination of the terms 

(06)2 = ^tV^z^-S^, 

The convergence of G n (x\, . . . ,x n ) to its continuum limit can now be studied in the 
effective theory, 

G n (x 1 ,... ,x n ) = ((j) (xi)...(j)o(x n )) con 

-a J d 4 y <0o(*i) . . . o (x n )£i(y)) con (1.2.13) 

n 

+a ^2 (M x i) ■ ■ ■ M x k) ■ ■ ■ M x n)) con + 0(a 2 ), 

k=l 

where the expectation values on the right-hand side are to be taken in the continuum 
theory with action Sq. 

Using the field equations 

For most applications, it is sufficient to compute on-shell quantities such as particle 
masses, S-matrix elements and correlation functions at space-time arguments, which are 
separated by a physical distance. It is then possible to make use of the field equations 
to reduce first the number of basis fields in the effective Lagrangian L\ and, in a second 
step, also in the 0(a) counter-term </>i of the effective composite fields. 

Of course, in the quantum theory, the field equations have to be used with care. 
Performing changes in the integration variables of the path integral shows immediately 
that the field equations are only valid up to additional terms. In particular, if one uses 
the field equations in the Lagrangian C\ under the space-time integral in eq. (1.2.13), 
the errors made are contact terms that arise when y comes close to one of the arguments 
Xi, . . . ,x n . Using the operator product expansion and the dimensions and symmetries, 
one easily verifies that these contact terms must have the same structure as the insertions 
of 4>i in the last term of eq. (1.2.13). The use of the field equations in C\ therefore just 
means a redefinition of the coefficients in the counter-term <pi (apart from contributions 



(1.2.12) 
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that are absorbed in the renormalization factor Z^, which originate from the mixing of 
(f>i with (p). 

It turns out that one may eliminate two of the terms in eq. (1.2.9) by using the field 
equations. A possible choice is to stay with the terms Oi, O3 and O5, which yields the 
effective continuum action for on-shell quantities to order a. Having made this choice 
one may apply the field equations once again to simplify the term 0i in the effective 
field as well. In the example of the axial current it is then possible to eliminate the 
term C 6 in eq. (1.2.12). 

1.2.2.4 Improved lattice action and fields 

The on-shell 0(a) improved lattice action is obtained by adding a counter-term to 
the unimproved lattice action such that the effects of the action S\ in the effective 
theory are canceled in on-shell amplitudes. This can be achieved by adding lattice 
representatives of the terms 0\, O3 and O5 to the unimproved lattice Lagrangian, with 
coefficients that are functions of the bare coupling g$ only. Leaving the discussion 
of suitable improvement conditions for later, we here note that the fields O3 and O5 
already appear in the unimproved theory and thus merely lead to a re-parametrization 
of the bare parameters go and uiq. In the following, we will not consider these terms 
any further. Their relevance in connection with massless renormalization schemes is 
discussed in detail in Ref. [19] . 

We choose the standard discretization F„ u of the field tensor [19] and add the 
improvement term to the Wilson-Dirac operator, eq. (1.1.43), 

Ampr = D + c sw —a^F^. (1.2.14) 

With a properly chosen coefficient c sw (<?o), this yields the on-shell 0(a) improved lattice 
action which has first been proposed by Sheikholeslami and Wohlert [42]. 

The perturbative expansion of c sw reads c sw = 1 + ciw<?o + 0(So)> with [43,44] 
Csw = 0.26590(7) where the Wilson plaquette gauge action is assumed. 

The 0(a) improved isospin currents and the axial density can be parametrized as 
follows, 

(4)J = ^a(1 + 6 A am q ){^ + ac A cVP a }, 

(Va)£ = Zy(l + b v am q ){V^ + ac v d u T^}, 

(P R ) a = Z P (l + b P am q )P a , (1.2.15) 

where T^=i^4r a V. 

We have included the normalization constants Za.v,p> which have to be fixed by ap- 
propriate normalization conditions. Again, the improvement coefficients 6a.v,p an d 
ca,v are functions of go only. At tree level of perturbation theory, they are given by 
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6a = bp = by = 1 and ca = cy = [19, 45] and to 1-loop accuracy and with the 
plaquette gauge action we have [43,46] 



c A = -0.005680(2) x^Cp + O^), 

c v = -0.01225(1) x glC ¥ + O(^), 

6 A = 1 + 0.11414(4) xg$C F + 0(g%), 

by = 1 + 0.11492(4) x^Cp + O^), 

b P = 1 + 0.11484(2) xglC ¥ + 0{gt). 



(1.2.16) 



Here Cp = 4/3. Non-perturbative determinations will be mentioned below. 
1.2.3 The PCAC relation 

We assume for the moment that on-shell O(a) improvement has been fully implemented, 
i.e. the improvement coefficients are assigned their correct values. If O denotes a renor- 
malized on-shell O(a) improved field localized in a region not containing x, we thus 
expect that the PCAC relation, 



be improved for this statement to be true. To see this we use again Symanzik's local 
effective theory and denote the O(a) correction term in c g- by (f>i. Eq. (1.2.17) then 
receives an order a contribution 



which is to be evaluated in the continuum theory. The PCAC relation holds exactly in 
the continuum and the extra term (1.2.19) thus vanishes, a conclusion that holds also 
in the Schrodinger functional. 

1.2.4 Non-perturbative improvement 

The coefficients of the different improvement terms need to be fixed by suitable improve- 
ment conditions. One considers pure lattice artifacts, i.e. combinations of observables 
that are known to vanish in the continuum limit of the theory. Improvement conditions 
require these lattice artifacts to vanish, thus defining the values of the improvement 
coefficients as a function of the lattice spacing. 

In perturbation theory, lattice artifacts can be obtained from any (renormalized) 
quantity by subtracting its value in the continuum limit. The improvement coefficients 
are unique and some of them have been cited above. 

Beyond perturbation theory, one wants to determine the improvement coefficients 
by MC calculations and it requires significant effort to take the continuum limit. It 





(1.2.19) 
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is therefore advantageous to use lattice artifacts that derive from a symmetry of the 
continuum field theory that is not respected by the lattice regularization. One may 
require rotational invariance of the static potential V(r), e.g. 



V(r = (2, 2, l)r/3) - V(r = (r, 0, 0)) = 



or Lorentz invariance 



[E(p)} 2 - [E(0)} 2 - p 2 = 



for the momentum dependence of a one-particle energy E. 

For 0(a) improvement of QCD it is advantageous to use violations of the PCAC 
relation (1.2.7), instead. PCAC can be used in the context of the Schrodinger functional 
(SF), where one has a large flexibility to choose appropriate improvement conditions 
and can compute the improvement coefficients also for small values of the bare coupling 
go, making contact with their perturbative expansions. A further - and maybe the 
most significant - advantage of the SF in this context is the following. As we have seen 
earlier, in the SF we may choose boundary conditions such that the induced background 
field, has non- vanishing components F^ v . Remembering eq. (1.2.14), we observe that 
correlation functions are then sensitive to the improvement coefficient c sw already at 
tree level of perturbation theory. With a vanishing background field this would be the 
case at higher orders only. Also in the non-perturbative regime this is the basis for a 
good sensitivity of the improvement conditions to c sw . 

As a consequence of the freedom to choose improvement conditions, the resulting 
values of improvement coefficients such as c sw , ca depend on the exact choices made. 
The corresponding variation of c sw , ca is of order a. It changes the effects of order a? in 
physical observables computed after improvement. We will come back to this point, but 
first we proceed to sketch the non-perturbative calculation of the coefficient c sw [19,43]. 

We define a bare current quark mass, m, viz. 



with the pseudo scalar boundary operator O a from eq. (1.1.39). When all improve- 
ment coefficients have their proper values, the renormalized quark mass, defined by the 
renormalized PCAC-relation, is related to m by 



It now suffices to choose 3 different versions of eq. (1.2.20) by different choices of the time 
coordinate xq and/or boundary conditions and obtain 3 different values of m, denoted 
by mi, m-2, m%. Since the prefactor in front of m in eq. (1.2.21) is just a numerical factor, 
we may conclude that all m^ have to be equal in the improved theory up to errors of 




(1.2.20) 




(1.2.21) 
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Figure 6: Unrenormalized current quark mass m in the improved theory, with non- 
perturbatively determined c sw and function of the time xq on a 32 x 16 3 lattice 

at f3 = 6.2 (a « 0.07fm)and k = 0.1350. The width of the corridor bounded by the 
dotted horizontal lines is 4 MeV. 

order a 2 . c sw and ca may therefore be computed by requiring 10 

nil = m 2 = m 3 ■ (1.2.22) 

This simple idea has first been used to compute c sw and ca as a function of go in the 
quenched approximation [19]. A good accuracy has been reached in these determinations 
for a < 0.1 fm. We comment on more recent determinations below, after showing some 
tests of the effectiveness of 0(a)-improvement, namely the test whether cutoff effects 
are indeed significantly reduced in practical calculations. 

Verification of O(a) -improvement 

The first test of the size of residual 0(a 2 ) effects is provided again by the PCAC 
relation. To set the scale, remember that the cutoff effects in the PCAC mass m were 
as large as several tens of MeV before improvement (cf. fig. 3). The situation after 
improvement, and for a somewhat larger value of the lattice spacing, is illustrated in 
fig. 6. Away from the boundaries, where the effect of states with energies of the order 
of the cutoff induces noticeable effects, m is independent of time to within ±2 MeV. 
Compared to the situation before O(a) improvement, this is a change by more than an 
order of magnitude. 

The second test is in the scaling of the vector meson mass, which is improved 
dramatically in the quenched approximation. In the range 0.01 < o? /r 2 . < 0.035 (a < 

10 In the practical calculations actually four different masses rrn were used to compute c sw and two 
more to extract ca- 
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0.1 fm) 11 and with statistical errors of 1-3%, no a effect was found (see Fig. 3, Ref. [33]). 
This is a significant test because i) the a-effects were large before improvement, ii) the 
statistical precision is good and iii) a 2 varies by more than a factor 3. 

In the quenched approximation several other tests of cutoff effects in the improved 
theory have been carried out [33, 48] revealing that 0(a)-improvement works. More 
precisely we mean by this that for a < 0.1 fm cutoff effects found after improvement are 
generally reasonably small (0(5%)) and well compatible with an 0(a 2 ) behavior. Still, 
no miracles have been achieved: at a ~ 0.1 fm also cutoff effects of about 15% have been 
found [33]. 

We now proceed to discuss a rather relevant detail in the non-perturbative deter- 
mination of improvement coefficients. 



1.2.4- 1 The constant physics condition 

We have mentioned before that beyond perturbation theory the improvement coefficients 
are not unique; they depend on the improvement condition. This ambiguity is of 0(a) 
in the coefficients and then of 0(a 2 ) in the physical observables after improvement. So 
everything is correct, but if these ambiguities are large, one has to take extra care. 

To give an explicit example, consider the quenched approximation but with non- 
degenerate quark masses, ^ m R,j- I n t ms theory the renormalized, improved 
currents and axial densities are given by a straight forward generalization of what we 
wrote down earlier [49], 

(Ar)J' = ZAil + bA^am^ + am^U^ + acAd^}, 

(P R )^ = Z P (1 + 6 P Karris + am^PV , (1.2.23) 

and the bare quark masses, moi, of the Lagrangian are related to the renormalized ones 
via 12 



i = — Z m q i(l + a 6 m m qii ) , m q i = m i - m c . (1.2.24) 



One then derives the (unrenormalized) PC AC relation [49] 
duAV + adlduP ij , 

— [1 + a 2 (6a - bp) (m q ,j + m qj )J 

= Z [m q ,j + m qJ + a b m (m 2 j + m 2 j)] + 0(a 2 ) . (1.2.25) 

This operator identity has to be understood as in eq. (1.2.17). Applying it with a few 
external operators/kinematical conditions, one can extract bA — bp, b m and Z [49,50]. 



11 We use the reference scale ro defined in terms of the static quark potential. [47] It has a phenomeno- 
logical value of ro ~ 0.5 fm. 

12 In the improved action, b m is the coefficient of the field O5. 
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Figure 7: The improvement coefficient b& — bp extracted on a 12 x 8 3 lattice and a lattice 
spacing of a ~ 0.1 fm. Axial current and density are inserted at distance xq. Squares 
are obtained with standard lattice representatives of the derivatives and 9*8^, while 
circles come from improved (and therefore less local) derivatives [51]. 



Suitable correlation functions have been implemented in the Schrodinger functional . 
As shown in Fig. 7, the result for b^ — bp at a ~ 0.1 fm does depend a lot on the details, 
namely here xq and the exact lattice representation of derivatives. This is an extreme 
example of an improvement coefficient that is difficult to determine. The reasons are 
surely that the 0(a)-effects are not dominating over the 0(a 2 ) effects in the considered 
correlation functions. Presumably the 0(a)-effects are just not very large. We proceed 
to discuss what one should do in these cases. 

When the improvement conditions are formulated through Schrodinger functional 
correlation functions, we can also study them in perturbation theory. One chooses 
kinematical variables such as T/L ,xq/L exactly as in the non-perturbative set up. One 
then computes the expansion of an improvement coefficient, here denoted generically 
by b, as a series 

b(g ,a/L) ~ &(°)(a/L) + g 2 b^(a/L) + O(<7 4 ) . (1.2.26) 

Taking the case of 6a — bp and L/a = 8,T/a = 12 as in Fig. 7, the perturbative coef- 
ficients show a similar dependence on the kinematics as the non-perturbative 
results [51]. 

While generically the functions b^°\a/L), b^\a/L) depend on the kinematical 
choices made in the improvement coefficients, the values b^°'(0), frW(0) are unique. 
This is the precise meaning of our earlier statement that improvement coefficients are 
unique in perturbation theory. 
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On the left we illustrates what 
happens generically. The contin- 
uous lines correspond to two dif- 
ferent improvement conditions at 
constant physics. By this we mean 
that all length scales in the correla- 
tion functions which define the im- 
provement condition are kept fixed 
in physical units, for example in 
units of ro- When go is varied, 
only the lattice spacing a changes 
and this is the situation to which 
the SET can be applied. The two 
different improvement coefficients 
then very rapidly go to an almost 
unique function of go as go — ► 0. 
On the other hand, one may want to set improvement conditions at fixed L/a for 
practical reasons 13 . In this setting, the 0(a) ambiguity in the improvement coefficients 
does not go to zero as go — > (and a — ► 0) (dotted lines). 

If the dependence of b^ l \a/L) on a/L is rather weak, say 

|6 (i) (a/L) -b«(0)| < 1(T\ (1.2.27) 

for the relevant a/L, one may also choose fixed a/L. Still it is clear from our discussion 
that fixed L/ro is to be preferred whenever possible. In cases where it is advantageous 
to work at finite quark masses, also the combination ro m^j should be kept constant. 
Note that these constant physics conditions do not have to be satisfied very precisely 
since we are talking about a correction to an 0(a) term. 

Such conditions have first been imposed in Ref. [51] - exactly because of Fig. 7. 
We show the result for 6a — bp, with the constant physics condition implemented, in 
Fig. 8. In this extreme case the ambiguity is of the same order as the improvement 
coefficient itself. In such a situation it is rather tempting to just put the coefficient to 
zero. However, it should be obvious that it is then not guaranteed that linear a-effects 
are absent after improvement. An extrapolation to the continuum using an 0(a 2 ) model 
function for the cutoff effects might then give significantly wrong results. We repeat this 
relevant fact in different words. While an unfortunately chosen improvement condition, 
but implemented at constant physics, may even enlarge the cutoff effects for intermediate 
a, it guarantees 14 that linear a-effects are absent. Putting the improvement coefficient 

13 In fact this was done originally for c sw , ca [52]. Note, however, that the small tree-level a-effect was 
subtracted from the non-perturbative ones to insure the improvement coefficients go to their tree-level 
values exactly and further the conditions were chosen such that also \b^\a/L) — b^(0)\ is very small. 

14 Here we use a strong wording. Remember of course that everything is based on the SET, which has 
not been proved. Remember also Sect. 1.2.2.2. 
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Generic behavior of improvement coefficients. 
Dotted lines correspond to two different im- 
provement conditions at fixed a/L, while the 
full lines are for fixed L in physical units. 
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Figure 8: Non-perturbative 6a — bp for iVf = and plaquette gauge action (left) as well 
as its ambiguity A (6a — bp). The latter is the difference of 6 a — bp for one improvement 
condition and 6a — bp for another one. 

to zero or some perturbative approximation does not guarantee the latter and the linear 
effects should at least be estimated in some way. 

Before we review what is known in 0(a)-improvement, we return once more to 
the improvement coefficients c sw ,ca- With plaquette gauge action, a/L = 1/8 and the 
chosen improvement condition, eq. (1.2.27) is very well satisfied for i = 0, 1 [43]. It 
appears justified to subtract these effects perturbatively as it was done in [52,53]. For 
ca there is evidence for a significant O(a) ambiguity [50,54,55]. As a consequence, 
the computation with Nf = 2 has been carried out with an improvement condition at 
constant physics [56]. 

1.2.4-2 Summary of results 

A summary of the available results for O(a) improvement coefficients is listed in Table 1. 
Many investigations have been carried out accumulating a quite advanced knowledge. 
It is impossible for me to review it all, rather I add only a few comments. 

• In the quenched approximation there is a rather complete knowledge for the Wil- 
son gauge action, but ca has not been determined with a constant physics condi- 
tion. Conservative continuum extrapolations in the improved theory usually start 
only with values of a which are a bit below 0.1 fm. 

The problem with "exceptional configurations" [64] , which is specific to the quenched 
approximation, is enhanced by c sw > [65]. Thus small quark masses can only be 
reached near the continuum limit. 

• For Nf = 2 and Wilson gauge action, the determination of improvement coeffi- 
cients is quite advanced. Interesting computations of physical observables can be 
done in the improved theory. 
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improvement coefficient 
or improved fields 


order 
in PT 




gauge action 


reference 




1-loop 




Wilson 
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1-loop 




Wilson 
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-ya Aa 
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1-loop 




Wilson 


[46, 58] 




1-loop 
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b m . 6a — bp, 6v, ca, cy 










t fa A a 
A* ' M 


JNP 
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Wilson 
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INF 
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Wilson 


[54] 
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1NP 
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Wilson 


[49] 


b A - b P ,b m 


NP at CP 





Wilson 


[51] 


6 A - o P ,6 m 


NP at CP 


t) 


Wilson 


[60] 


Csw 


NP 
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Wilson 


[53] 


Csw 


NP 
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Wilson 


[61] 


CA 


NP at CP 
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Wilson 


[56] 


b m , bA, bp, by, Cy, b T 










ya j^a rpa 


1-loop 




1 x 1 & 1 x 2 


[58] 


Csw 


1-loop 




1 x 1 & 1 x 2 


[62] 


Csw 


NP at CP 


0,2,3 


Iwasaki 


[63] 


CA 


NP at CP 


3 


Iwasaki 


CP-PACS & JLQCD 



Table l: Literature on 0(a) improvement coefficients. NP stands for non-perturbative, 
CP for constant physics condition, Wilson for the Wilson plaquette action, "1 x 1 & 
1 x 2" for the gauge action including a 1 x 2 loop and Iwasaki for the action where the 
coefficient of that term is set to the value proposed by Y. Iwasaki. 

• For iVf = 3 and Wilson gauge action, evidence for the existence of a first order 
phase transition in the K)-plane has been reported [66]. The value of c sw was 
fixed to typical values found for iVf = 0,2. The authors conclude that one has 
to remain with lattice spacings significantly below a = 0.1 fm for this action. 
Recent investigations [67] have shown that simulations with Wilson fermions at 
small quark masses are only algorithmically stable when the physical volume is 
sufficiently large. 15 It is possible that the simulations that lead to the evidence 
for a phase transition [66] suffer from this problem. We therefore consider the 
existence (and even more its position) of this phase transition not as settled. 
Unfortunately it appears that only large volume simulations with good statistics 
can tell. Note that this statement may be controversial. 

• With Iwasaki gauge action, c sw is now known [63] non-perturbatively for iVf = 
15 The lattice spacing should also be sufficiently small, but this is less relevant here. See also [68]. 
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0, 2, 3. Non-perturbative computations of ca are in progress and the relevant 
other improvement coefficients are known to 1-loop order [58]. A small grain of 
salt is the following. The perturbative expansion seems less trustworthy for this 
action, since the bare coupling is significantly larger than for the Wilson gauge 
action at the same lattice spacing, while the perturbative coefficients are of a 
similar size. 

• In Sect. 1.2.4.4 we will come to the question whether more information is needed 
to do interesting computations in QCD. 

1.2. 4- 3 0(a 2 ) effects after improvement for 2Vf = 2 

The ALPHA collaboration has performed a large series of simulations with the Schrodinger 
functional in small volumes of a linear extent of at most L ~ 0.5 fm [23,69,70]. The 
purpose of these simulations was the non-perturbative determination of the running of 
scale dependent quantities (see Sect. II. 1). Extracting the continuum limit of the step 
scaling functions would have been impossible without 0(a)-improvement. It is fair to 
say that 0(a)-improvement worked marvelously well in this situation. The interested 
reader may in particular study the discussion in [69]. The emerging picture for QCD 
in the Wilson formulation is that the Symanzik expansion and improvement works very 
well for lattice spacings around 0.05 fm and below. 

At larger lattice spacings, the situation is unclear at present. In particular indirect 
evidence has been shown that 0(a 2 ) effects are rather large at a ~ 0.1 fm [22,71]. 
Large scale simulations in large volume will tell, whether the remaining 0(a 2 ) effects 
are a problem in practice. We emphasize that there is no evidence anywhere that the 
Symanzik theory of o-effects is invalid. The remaining issue is how small a has to be 
made for it to be precise when truncated to the first one or two correction terms. It is 
also expected that smaller lattice spacings are needed for smaller quark masses. This 
issue is discussed in Steve Sharpe's lectures. [72] 

1.2.4-4 Do we need more? 

The question arises, whether the knowledge described above is sufficient to do QCD 
computations. 

First, we recall that we assumed that one is interested in on-shell quantities. In 
the popular MOM-scheme [73] used for the renormalization of composite operators, one 
does, however, consider off-shell correlation functions. While improvement in that situ- 
ation has been investigated, too [74-76], we continue under the assumption that off-shell 
correlation functions are avoided. So far we have restricted our discussion to the case 
of QCD with mass-degenerate quarks. This does of course not correspond to real QCD. 
General quark masses are no complication for quenched quarks; we already referred 
to that situation in Sect. 1.2.4.1. However, dynamical quarks of different masses allow 
for a whole series of new O(a) counter-terms which are proportional to the difference 
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of quark masses [77]. Determining them non-perturbatively would be a considerable 
challenge. However, let us look at their influence in practice 16 . For TVf = 3, the biggest 
mass difference is of order m s . For a < 0.1 fm, we have 

am^ s < 1/30. (1.2.28) 

Thus we are talking about very small effects even at the largest typical lattice spacing. 
It is then sufficient to know the coefficient with, an accuracy of 0.25 or so which is surely 
possible by perturbation theory. The same argument applies to b g [19]. Note, however, 
that for non-degenerate masses, eq. (1.2.24) is modified already at the level a and this 
has to be taken into account [77], or better it is avoided by relying on the PCAC masses 
only. 

Obviously, for a dynamical charm quark, the question should be reconsidered care- 

fully. 

Another issue is the improvement of more complicated composite operators, such 
as 4-fermion operators. Many improvement terms may be necessary. For these difficult 
cases one may consider a mixed action approach, where (some part) of chiral symmetry 
forbids the 0(a) terms [78]. Also twisted mass lattice QCD offers an interesting strategy 
for these problems. We refer to Stefan Sint's lectures. 

In conclusion, a very large class of interesting problems do not seem to need more 
beyond what is presently known or is being determined. One may therefore take ad- 
vantage of the simplicity of the formulation and work on the algorithms [79-84] to be 
able to work at small lattice spacings and quark masses; see also the lectures of Tony 
Kennedy [85]. 



We thank M. Liischer for emphasizing this point. 
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Lecture II. 

Fundamental parameters of QCD from the lattice 



II. 1 The problem of scale dependent renormalization 

Let us investigate the extraction of short distance parameters (Sect. 1.2a) in some detail. 
After some brief comments on the conventional way of obtaining a from experiments 
we explain how one can compute it at large energy scales using lattice QCD. This paves 
the road for a computation of the A-parameter and then also for the renormalization 
group invariant quark masses. We will finally comment briefly on other scale-dependent 
renormalization. 



II. 1.1 The extraction of a from experiments 

One considers experimental observables depending on an overall energy scale q and 
possibly some additional kinematical variables denoted by y. The observables can be 
computed in a perturbative series which is usually written in terms of the MS coupling 
a— 17 

"MS' 

Oi(q,y) = a m (n) + Ai(y)a^{n) + . . . , H = q. (11.1. 1) 

For example Oi may be constructed from jet cross sections and y may be related to the 
details of the definition of a jet. 

The renormalization group describes the energy dependence of a in a general scheme 

^ = m > (n.i.2) 

where the /3-function has an asymptotic expansion 

(3(g) -g 3 {b + g 2 h + ...} , (II.1.3) 

b = (^(H - |iV f ) , h = ^(102 - f N { ) , 

with higher order coefficients bi, i > 1 that depend on the scheme. Eq. (II. 1.3) entails 
asymptotic freedom: at energies that are high enough for eq. (II. 1.3) to be applicable 
and for a number of quark flavors, iVf , that is not too large, a decreases with increas- 
ing energy as indicated in Fig. 9. The solution of eq. (II. 1.2) contains an integration 
constant, the renormalization group invariant parameter A. It is (exactly) given by 



17 We can always arrange the definition of the observables such that they start with a term a. For 
simplicity we neglect all quark mass dependences; they are irrelevant for the main points of the present 
discussion. 
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Figure 9: Running of g 2 and its reconstruction from the step scaling function (points). 

where g = g(fJ.) Note that A is different in each scheme. If a coupling ax is related to 
another one ay at the same energy scale via 

ay (At) =ax{fi) + cf Y [a x (fi)} 2 + c^ Y [a x (/i)] 3 + • ■ ■ , (H.1.5) 

their A-parameters are converted via 

Ax/ Ay = exp{-cf y /(8^ )} • (H.1.6) 

From the above equations it is easy to show that eq. (II. 1.6) is exact. For large \i one 
reads off the asymptotics 

-2^00 1 MnH/^/A 2 )] / {ln[ln( M 2 /A 2 )]} 2 \ 

9 ~ 6 ln(/i 2 /A 2 ) fc3[ln(/i 2 /A 2 )] 2 + \ [ln(^ 2 /A 2 )] 3 )' K ' ' 

We note that - neglecting experimental uncertainties - extracted in this way 
is obtained with a precision given by the terms that are left out in eq. (II. 1.1). In 
addition to a 3 -terms, there are non-perturbative contributions which may originate from 
"renormalons" , "condensates" (the two possibly being related), "instantons" or - most 
importantly - may have an origin that no physicist has yet uncovered. Empirically, one 
observes that values of determined at different energies and evolved to a common 
reference point using the renormalization group equation eq. (II. 1.2) including 62 agree 
rather well with each other; the aforementioned uncertainties are apparently not very 
large (Fig. 10). Nevertheless, determinations of a are limited in precision because of 
these uncertainties and in particular if there was a significant discrepancy between a 
determined at different energies it would be difficult to say whether this was due to 
the terms left out in eq. (II. 1.1) or was due to terms missing in the Standard Model 
Lagrangian. 

It is an obvious possibility and at the same time a challenge for lattice QCD to 
achieve a determination of a in one (non-perturbatively) well defined scheme and evolve 
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Figure 10: The running coupling in the MS scheme extracted from various scattering 
experiments compared to the perturbative scale dependence. Graph of the particle data 
group [86]. 



this coupling to high energies. There one may use eq. (II. 1.4) with a perturbative 
approximation for /3(g). For a good precision 62 should be known. 

The A-parameter can then serve as an input for perturbative predictions of jet cross 
sections or the hadronic width of the Z-boson and compare to high energy experiments 
to test the agreement between theory and experiment. Since in the lattice regulariza- 
tion, QCD is naturally renormalized through the hadron spectrum, such a calculation 
provides the connection between low energies and high energies, verifying that one and 
the same theory describes both the hadron spectrum and the properties of jets. 

Note. A dis-satisfying property of a^g is that it is only defined in a perturbative 
framework; strictly speaking there is no meaning of phrases like "non-perturbative cor- 
rections" in the extraction of a^jg from experiments. The way that I have written 
eq. (II. 1.1) suggests immediately what should be done instead. An observable Oi itself 
may be taken as a definition of a - of course with due care. Such schemes called physical 
schemes are defined without ambiguities. This is what will be done below for observ- 
ables that are easily handled in MC-simulations of QCD. For an additional example 
see [87,88]. 

II. 1.2 Reaching large scales in lattice QCD 

Let us simplify the discussion and restrict ourselves to the pure Yang- Mills theory with- 
out matter fields in this section. A natural candidate for a non-perturbative definition 
of a is the following. Consider a quark and an anti-quark separated by a distance r and 
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in the limit of infinite mass. They feel a force F(r), the derivative of the static potential 
V(r), which can be computed from Wilson loops (see e.g. [89]). A physical coupling is 
defined as 

MA*) = 7^r 2 F(r), fx = 1/r, C F = 4/3 . (II.1.8) 

It is related to the MS coupling by eq. (II. 1.5) with a certain constant c^ Sq ^, which also 
determines the ratio of the A-parameters (eq. (II. 1.6)). Note that G^qq IS cL renormalized 
coupling defined in continuum QCD. 

Problem. If we want to achieve what was proposed in the previous subsection, the 
following criteria must be met. 

• Compute ce q q(/i) at energy scales of /i ~ lOGeV or higher in order to be able to 
make the connection to other schemes with controlled perturbative errors. 

• Keep the energy scale \x removed from the cutoff a" 1 to avoid large discretization 
effects and to be able to extrapolate to the continuum limit. 

• Of course, only a finite system can be simulated by MC. To avoid finite size effects 
one must keep the box size L large compared to both the mass of the lightest 
physical state (the pion) as well as to a typical QCD scale, say the potential scale 
r [47]. 

These conditions are summarized by 

L » r , — ~ 1 » - ~ — ?T ^ 7 » a , (II.1.9) 
TOy,- 0.14CeV fj, lOGeV 

which means that one must perform a MC-computation of an A 4 lattice with N = 
L/a 3> 70. In the near future it is impossible to perform such a computation. The origin 
of this problem is simply that the extraction of short distance parameters requires that 
one has control over physical scales that are quite disparate. To cover these scales in one 
simulation requires a very fine resolution, which is too demanding for a MC-calculation. 

Of course, one may attempt to compromise in various ways. E.g. one may perform 
phenomenological corrections for lattice artifacts, keep l//i ~ a and at the same time 
reduce the value of \i compared to what I quoted in eq. (II. 1.9). Calculations of a qq - 
along these lines have been performed in the Yang-Mills theory [90-92]. It is difficult to 
estimate the uncertainties due to the approximations that are necessary in this approach. 
More recently, results in the continuum limit could be obtained up to /i = 4 GeV by 
simulating very large lattices [93,94], and still it is not obvious that one has reached the 
perturbative region. We will come back to this. 
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Solution. Fortunately these compromises can be avoided altogether [37]. The solution 
is to identify the two physical scales, above, 

H = 1/L . (II.1.10) 

In other words, one takes a finite size effect as the physical observable. The evolution 
of the coupling with [i can then be computed in several steps, changing /i by factors of 
order 2 in each step. In this way, no large scale ratios appear and discretization errors 
are small for L/a S> 1. 

The complete strategy to compute short distance parameters is summarized in 
Fig. 11. One first renormalizes QCD replacing the bare parameters by hadronic observ- 
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Figure ll: The strategy for a non-perturbative computation of short distance parame- 
ters. SF refers to the Schrodinger functional scheme introduced in Sect. 1.1.6 

ables. This defines the hadronic scheme (HS) as explained in Sect. 1.1. At a low energy 
scale fj, = 1/ L max this scheme can be related to the finite volume scheme denoted by SF 
in the graph. Within this scheme one then computes the scale evolution up to a desired 
energy fi = 2 n /L max . As we will see it is no problem to choose the number of steps n 
large enough to be sure that one is in the perturbative regime. There perturbation the- 
ory (PT) is used to evolve further to infinite energy and compute the A-parameter and 
the renormalization group invariant quark masses. Inserted into perturbative expres- 
sions these provide predictions for jet cross sections or other high energy observables. 
In the graph all arrows correspond to relations in the continuum; the whole strategy 
is designed such that lattice calculations for these relations can be extrapolated to the 
continuum limit. 

For the practical success of the approach, the finite volume coupling (as well as the 
corresponding quark mass) must satisfy a number of criteria. 
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• They should have an easy perturbative expansion, such that the /3-function (and 
r-function, which describes the evolution of the running masses) can be computed 
to sufficient order. 

• They should be easy to calculate in MC (small variance!). 

• Discretization errors must be small to allow for safe extrapolations to the contin- 
uum limit. 

Careful consideration of the above points led to the introduction of renormalized cou- 
pling and quark mass through the Schrodinger functional (SF) of QCD [12,17,18,32,95], 
introduced in Sect. 1.1. In the Yang-Mills theory, an alternative finite volume coupling 
was studied in detail in [96,97]. 

The criteria eq. (II. 1.9) apply quite generally to any scale dependent renormaliza- 
tion, e.g. the one of 4-fermion operators of the effective weak Hamiltonian at scales 
fi <C M\y- Indeed, details have been worked out for several cases [9-11,98-101]. 

A frequently applied alternative is to search for a "window" where f/, is high enough 
to apply PT but not too close to a" 1 , the cutoff [102]. An essential advantage of 
the details of the approach of Ref. [102] as applied to the renormalization of compos- 
ite quark operators is its simplicity: formulating the renormalization conditions in a 
MOM-scheme 18 , one may use results from perturbation theory in infinite volume in the 
perturbative part of the matching. Since, however, it is difficult to reach high energies 
fi in this approach, we will not discuss it further and refer to [102-105] for an idea of 
the present status and further references, instead. 

II. 2 The computation of a((j,) and A 

We are now in the position to explain the details of Fig. 11 [37,98]. The problem has been 
solved "completely" in the SU(iV) Yang-Mills theories for N = 2,3 [17,93,95,97,98,106]. 
In QCD with only two dynamical quarks the strategy has been carried out well [69,70], 
except for the last line in the graph, which needs more work. In the present context, 
the Yang-Mills theory is of course equivalent to the quenched approximation of QCD 
or the limit of zero flavors. We will therefore also refer to results in quenched QCD. 

Our central observable is the step scaling function that describes the scale-evolution 
of the coupling, i.e. moving vertically in Fig. 11. The analogous function for the running 
quark mass will be discussed in the following section. 

II. 2.1 The step scaling function 

We start from a given value of the coupling, u = g 2 {L). When we change the length 
scale by a factor s, the coupling has a value g 2 (sL) = u'. 19 The step scaling function, 

18 We do not use the prefix RI (regularization independent), since also the SF-scheme is regularization 
independent. 

19 Figure 11 is for the most frequently used case s = 2. 
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Figure 12: The computation of a lattice step scaling function. 



a is then defined as 



a(s, u) 



i 

u . 



(II.2.1) 



The interpretation is obvious. a(s, u) is a discrete /3-function. Its knowledge allows for 
the recursive construction of the running coupling at discrete values of the length scale, 



u k =g 2 (s k L) , 



(II.2.2) 



once a starting value uq = g 2 (L) is specified (cf. the points in Fig. 9). The step 
scaling function, a, which is readily expressed as an integral of the /^-function, has a 
perturbative expansion 



a(s,u) = u + 2&o ln(s)u 2 + . . . . 



(II.2.3) 



On a lattice with finite spacing, a, the step scaling function will have an additional 
dependence on the resolution a/L. We define 

E(s,u,a/L) = u' , (II.2.4) 

with 

g 2 (L) = u, g 2 {sL) = u , go fixed, L/a fixed . (II. 2. 5) 
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The continuum limit a(s,u) = T,(s,u, 0) is then reached by performing calculations for 
several different resolutions and extrapolation a/L — ► 0. The computation of a(2,u) is 
illustrated in Fig. 12. In detail, one performs the following steps: 

1. Choose a lattice with L/a points in each direction. 

2. Tune the bare coupling go such that the renormalized coupling g 2 (L) has the value 
u and tune the bare quark mass, too, such that the PCAC-mass, defined at fixed 
physical kinematics [69], vanishes. 20 

3. At the same value of go, simulate a lattice with twice the linear size; compute 
v! = g 2 (2L). This determines the lattice step scaling function £(2,u, a/L). 

4. Repeat steps 1.-3. with different resolutions L/a and extrapolate a/L — > 0. 

Note that step 2. takes care of the renormalization and 3. determines the evolution of 
the renormalized coupling. 

The presently most advanced numerical results are displayed in Fig. 13. The cou- 
pling used is exactly the one defined in Sect. 1.1.6 and the calculation is done in the 
theory with iVf = 2 flavors of 0(a)-improved fermions. One observes that the depen- 
dence on the resolution is very weak, in fact it is not observable within the precision of 
the data in Fig. 13. We now investigate in more detail how the continuum limit of S is 
reached. As a first step, we turn to perturbation theory. 



II. 2. 2 Lattice spacing effects in perturbation theory 

Symanzik has investigated the cutoff dependence of field theories in perturbation the- 
ory [107]. Generalizing his discussion to the present case, one concludes that the lattice 
spacing effects have the expansion 



£(2,u,a/L) - a(2,u) 



5 1 (a/L)u + 5 2 (a/L)u 2 + ... (II.2.6) 



a(2,u) 

n 

5 n (a/L) a/L ~° ^ efcin [ln(f)] fe (f)+4,n[ln(f)] fc (f) 2 + --- ■ 

We expect that the continuum limit is reached with corrections 0(a/L) also beyond 
perturbation theory. In this context 0(a/L) summarizes terms that contain at least 
one power of a/L and may be modified by logarithmic corrections as it is the case 
in eq. (II. 2. 6). To motivate this expectation recall Sect. 1.3, where we explained that 
lattice artifacts correspond to irrelevant operators which carry explicit factors of the 
lattice spacing. Of course, an additional a-dependence comes from their anomalous 
dimension, but in an asymptotically free theory such as QCD, this just is a logarithmic 
(in a) modification. 



The details of the definition of the kinematics affects only the cutoff effects, see Sect. 1.2.2. 
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Figure 13: The lattice step scaling function after 2-loop observable improvement for 
N f = 2. 
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Figure 14: Lattice artifacts at 1-loop order. Points show 5\{a/L) for the SU(3) Yang- 
Mills theory with 1-loop improvement. Wilson (circles), Symanzik (squares), Iwasaki 
(filled triangles) and DBW2 (empty triangles) are shown. The dotted line corresponds 
to the linear piece in a, when only tree-level improvement is used, instead (only for 
Wilson). 

As mentioned in Sect. 1.1.8, the lattice artifacts may be reduced to 0((a/L) 2 ) 
by canceling the leading irrelevant operators. In the case at hand, this is achieved 
by a proper choice of c t (<?o)- It is interesting to note, that by using the perturbative 
approximation 

c t (g ) = l + ci 1) g 2 (11.2.7) 

one does not only eliminate e^^n f° r n = 0,1 but also the logarithmic terms generated 
at higher orders (n > l)are reduced, 

e n ,n = 0, e n _i jn = . (II.2.8) 

For tree-level improvement, c t (<?o) = 1, the corresponding statement is e njn = 0. Heuris- 
tically, the latter is easy to understand. Tree-level improvement means that the prop- 
agators and vertexes agree with the continuum ones up to corrections of order 0(a 2 ). 
Terms proportional to a can then arise only through a linear divergence of the loop 
integrals. Once this happens, one cannot have the maximum number of logarithmic 
divergences any more; consequently e n ,n vanishes. 

To demonstrate further that the abelian field introduced in the previous section 
induces small lattice artifacts, we show 5\(a/L) for the one loop improved case in Fig. 14 
(circles). The term that is canceled by the proper choice = —0.089 is shown as a 
dashed line. The left over 0((a/L) 2 )-terms are below the 1% level for couplings u < 2 
and lattice sizes L/a > 6. For not too large coefficient of the 1x2 loop in the action, 
they are close to quadratic in a/L in the range of interest. Also the fermion contribution 
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Figure 15: The lattice step scaling function after 2-loop observable improvement for 
iVf = 0. Data from [98, 113]. The point at o = is the continuum limit of [98] taken 
with a subset of the data. 

and the 2-loop cutoff effects are known [16,57]. When the 0(a)-improved theory is used 
(and the coupling has reasonable values), they turn out to be smaller than the 1-loop 
terms discussed here. 

Of course the cutoff effects also depend on the gauge action. In one class of such 
actions, one adds a 1 x 2 rectangular loop. Its coefficient has originally been deter- 
mined to have tree-level Symanzik 0(a 2 ) improvement in the pure gauge theory [108]. 
Later other choices have been proposed by Y. Iwasaki [109] and the QCD-TARO col- 
laboration [110], based on renormalization group considerations. We refer to these as 
"Symanzik, Iwasaki, DBW2" . The choice of the action close to the Schrodinger func- 
tional boundaries does of course involve an additional freedom discussed in [111,112]. 
For the choice favored in [112], we included the corresponding 8\{a/L) in Fig. 14. For 
the tree-level Symanzik improved theory the cutoff effects are remarkably small. One 
might think that this is so by construction, but note that here we are discussing 1-loop 
effects and in addition also boundary operators contribute to the SET expansion of the 
cutoff effects. 

Altogether, we now understand better why the a/L-dependence is so small in 
Fig. 13. Still, because the continuum limit is so important, and we have seen that 
unexpected difficulties may be present (Sect. 1.2.2.2), we continue its discussion in the 
pure gauge theory. There, numerical simulations with better resolutions have been 
carried out. 

II. 2. 3 The continuum limit — universality 

In fact, it is not only of interest to investigate how exactly the continuum limit is 
approached. Its very existence and its universality, i.e. the independence of the renor- 
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Figure 16: Universality test in the SU(3) Yang Mills theory. The data from top (trian- 
gles) to bottom (open circles) are for the Iwasaki, Symanzik and Wilson gauge action. 
Both the boundary improvement of the action and the improvement of the observables 
are carried out. At present this is possible at the 2-loop level for the Wilson gauge 
action (same as Fig. 15) and at the 1-loop level otherwise. Data are from [69,114]. 



malized continuum observable from the lattice action, is worth testing by precise MC 
computation. 

Above we have discussed how one can add boundary terms to the action to system- 
atically reduce the cutoff effects linear in a. Given the knowledge of the observables in 
perturbation theory at finite lattice spacing, we can go further and remove the cutoff- 
effects for arbitrary a/L, i.e. for instance including 0((a/L) 2 ), but up to a fixed order 
in u. This is achieved by the improvement of the observable [69,70,97] 

eW(2, M /I) = ^1^13 , (II .2.9) 

The approach has been carried out with k = 2 for the Af = 2 data (Fig. 13) and is 
tested thoroughly for Af = in Fig. 15. 

Recently, a nice universality test has been carried out by the CP-PACS collabo- 
ration [114]. Obviously the results in Fig. 16 are well compatible with an a — > limit 
which is independent of the action. Similar results exist for the SU(2) theory [97] . These 
results leave little doubt that the continuum limit of the Schrodinger functional exists 
and is independent of the lattice action. In turn this also supports the statement that 
the Schrodinger functional is renormalized after the renormalization of the coupling 
constant. 

We return to the extraction of a in Af = 2 QCD. Even though a-effects are not 
statistically significant, one has to be careful how one extracts the continuum limit. The 
worry does not so much concern the central values, but the correct estimate of their 
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Figure 17: The QCD /^-function in the SF scheme. 



uncertainties. For example just averaging data at all values of a/L produces an unreal- 
istically small statistical error, because one has then assumed that a-effects are entirely 
absent, although the data tell only that they are smaller than statistical uncertainties. 
One possible strategy for the continuum extrapolation is thus a fit to a constant that 
uses the lattices with L/a = 6, 8 only. As a check of this procedure, different variants of 
a combined continuum extrapolation of all the data sets, but excluding L/a = 4 were 
carried out. For example the ansatz 



with a constant p was fitted to the data. The final conclusion was that the simple fit to 
a constant (for L/a = 6, 8) yields realistic error estimates for the existing data set [70]. 
A further check of this procedure is Fig. 15 where the dotted lines represent the error 
band obtained in this way and the four points at smaller lattice spacing are in perfect 
agreement with this band. 

After this long - but important - discussion of cutoff effects, we are convinced 
that we have continuum results for the step scaling function with realistic uncertainties. 
They are ready to be used to construct the running coupling and the A-parameter. 

II. 2. 4 The running of the coupling 

The numerical values of cr(u) are next represented by a smooth interpolating function 
(a polynomial in u). With this function the running coupling g 2 (2~ l L max ) = m can be 
constructed from the recursion 



£ (2) (2,u,a/L) = a(2,u)+pu 4 (a/L) 2 





= 0...n; 



(II.2.10) 
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the result is shown in Fig. 18 for the (arbitrary) choice 
a recursion for the (5- function itself [70], 



0{ 



^U i + l /Ui <j'(Ui + l) P(y/U i+1 ] 



5.5. One can also set up 



(II.2.11) 



Together with a start value for the /3-function taken from perturbation theory (3-loop) 
at the weakest coupling (a ~ 0.08) this yields the numerical results Fig. 17. Their 
agreement with perturbation theory is excellent at weak couplings a < 0.2, while at the 
largest couplings significant deviations from perturbation theory are present for Nf = 2. 
Indeed the difference between non-perturbative points and 3-loop can't be described by 
an effective 4-loop term with a reasonable coefficient. At the same time the perturbative 
series just by itself does not show signs of its failure at, say, a ~ 0.3: instead successive 
orders yield smaller and smaller corrections. 

We return to the running couplings shown in Fig. 18. In the zero flavor case, 
also the region of [i of around 250 MeV was investigated with a specifically adapted 
strategy [113]. In this region, the SF coupling shows the rapid growth expected from a 
strong coupling expansion. 

Initially, the graphs Fig. 18 are obtained for /j, in units of /U m i n = 1/L max . One 
chooses ii max relatively large, but within the range covered by the non-perturbative 
computation of o~(u). The artificial scale L max has been replaced by the A parameter 
by use of eq. (II. 1.4). We proceed to explain this step. 



i — i i 1 1 1 1 ii 1 — i i 1 1 1 1 ii 1 — i 1 1 1 1 1 



i i 1 1 1 1 1| 1 — i i 1 1 1 i 1 — i i 1 1 1 1 1| 

SF scheme, N f =2 - 
3-loop 



SF scheme, N f =0 





10 2 fj./A 10 a 



Figure 18: Running coupling for Nf = 2 (left) and Nf = (right). 
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11. 2. 5 The A parameter 

We may evaluate eq. (II. 1.4) for the last few data points in Fig. 18 using the 3- loop 
approximation to the /3-function in the SF-scheme. The resulting A-values are essentially 
independent of the starting point, since the data follow the perturbative running very 
accurately at large \x. This excludes a sizable contribution to the /3-function beyond 
3-loops in this region. Indeed, a typical estimate of a 4-loop term in the /^-function 
would change the value of A by a tiny amount. The corresponding uncertainty can be 
neglected compared to the statistical errors. 

Changing A from the SF scheme to the MS scheme we then have 

N{ = , n max = 3.48 : ln(A m L max ) = -0.84(8) , (II.2.12) 
N { = 2, u max = 4.61 : ln(A m L max ) = -0.40(7) . (II.2.13) 

It remains to connect the artificially defined length scale L max to an experimentally 
measurable low energy scale of QCD such as the proton mass or the kaon decay constant, 
F K - 

So far it has been convenient to first evaluate L max in units of the low energy scale 
tq. This reference scale is precisely defined through the QCD static quark potential [47] 
but related to experiments only through potential models: tq ~ 0.5 fm. For Nf = a 
detailed investigation resulted in the continuum limit [93, 106] 

Nf = , u max = 3.48 : L max /r = 0.738(16) (II.2.14) 

and thus 

A^r = 0.60(5) . (II.2.15) 

In the iVf = 2 theory, the situation is illustrated in Table 2 which relies on results 
for ro/a from [115, 116]. On the one hand, all the numbers in italic are consistent, 
indicating that lattice spacing effects are small, on the other hand the first column 
shows that a is not yet varied very much. At the moment 

Ag.ro = 0.62(4) (4) 

is quoted, where the second error generously covers the range of numbers in italic in 
Table 2 and the first one comes from eq. (II. 2. 13). 

11. 2. 6 Discussion 

The scale dependence of the SF coupling is close to perturbative below «sf = 0.2 
and becomes non-perturbative above osf — 0.25. In fact a strong coupling expansion 
suggests that this coupling grows exponentially for large L. In the Nf = theory it 
was possible to verify this behavior for L close to lfm [113] (Fig. 18). Apart from 
the determination of A, an achievement of this investigation is the confirmation that 
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^max 


= 3.65 


^max 


= 4.61 


P = 6/<7 2 


ro/a 


■^max/ ^ 


A MS r O 


^max/ & 


A MS r O 


5.00 
5.29 
5.40 


5.45(5)(20) 
6.01(4) (22) 
7.01(5)(15) 


4.00(6) 
4.67(6) 
5.43(9) 


0.655(27) 
0.619(25) 
0.621(17) 


6.00(8) 
6.57(6) 
7.73(10) 


0.610(25) 
0.614(24) 
0.609(16) 



Table 2: A-parameter in units of ro for Nf = 2 for different resolutions in the low energy 
part of the calculation. Two different values of u max are considered, but A^jg-ro should 
be independent of these when cutoff effects are small. 

the transition between this non-perturbative region and the perturbative one is very 
smooth. 

With some effort this transition region was also bridged for a qq - in the SU(3) pure 
gauge theory. As shown in Fig. 19, the non-perturbative continuum results agree reason- 
ably well with the perturbative prediction at the largest (i (small r). The prediction is 
parameter free, since the A-parameter, eq. (II. 2. 15), can be changed to the qq-scheme, 
eq. (II. 1.8), where also the 3-loop /3-function is known [94, 117, 118]. The figure also 
illustrates that a reliable determination of A is not possible from a q q because for the 
accessible \i the 3-loop contribution to the j3- function is still very significant. 



source 


N { =0 


A MS r 

Nf =2 Nf =4 


N { =5 


ALPHA [17,70] 


0.60(5) 


0.62(6) 




world average [119] 




0.74(10) 


0.54(8) 


DISNNLO [120] 




0.57(8) 





Table 3: Results for A^ro for different number of flavors. 

In Table 3 we compare our results for A^ro to selected phenomenological ones, 
where we set ro = 0.5 fm. There appears to be an irregular A^-dependence, but we note 
that 

1. the errors are not very small yet, 

2. the 4-flavor A is obtained from the 5-flavor one by perturbation theory [121]. For 
this to be accurate, perturbation theory has to be accurate for /j, <C mb ea uty> which is 
not completely obvious. 

II. 2. 7 Improvements are necessary 

The most urgently needed improvement of the present results is to eliminate the model 
dependence which is intrinsic in the use of ro- So one should replace L max /ro by -Fk x 
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Figure 19: Running coupling in the qq-scheme (eq. (II. 1.8)). The perturbative prediction 
is given by the relation eq. (II. 1.4) between the running coupling and a£= . The dotted 
lines show the uncertainty of A eq. (II. 2. 15). 



L mm (computed at small enough light quark masses and small a) and an effort is 
presently being made. Also the strange quark sea has to be included ("2+1") and one 
should estimate the effects of the charm quark. Such 2+1 simulations are for example 
being carried out by JLQCD and CP-PACS, who have also studied the computation of 
the SF coupling with the gauge actions they are using [114]. 

These improvements will come and I am convinced that lattice results will yield 
the best controlled and most precise results for A in the long run. The reason is simple: 
I essentially described all the sources of systematic errors. The kind of assumptions one 
has to make are minimal. 

Let me also mention that there is a large number of other results for a from lattice 
gauge theory, where a([i ~ a -1 ) is extracted from quantities related to the cutoff. Some 
of them cite a very small error [122]. As discussed in [31], it is not easy to estimate the 
systematic errors in these computations, mainly because one cannot separately discuss 
higher order perturbative corrections and discretization errors. We thus think it is very 
desirable to carry out the program of Fig. 11 with good precision and the relevant 
number of flavors. 
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II. 3 Renormalization group invariant quark masses 



The computation of running quark masses and the renormalization group invariant 
(RGI) quark mass proceeds in analogy to the computation of a(fx). Since we are using 
a mass-independent renormalization scheme (cf. Sect. 1.1.7.2), the renormalization (and 
thus the scale dependence) is independent of the flavor of the quark. When we consider 
"the" running mass below, any one flavor can be envisaged; the scale dependence is the 
same for all of them. 

The renormalization group equation for the coupling eq. (II. 1.2) is now accompanied 
by one describing the scale dependence of the mass, 

f)vn 

^ = r(g) , (II.3.1) 
where r has an asymptotic expansion 

r(g) -f {d + g 2 d 1 + ...}, d = 8/(4vr) 2 , (II.3.2) 

with higher order coefficients di, i > which depend on the scheme. 

Similarly to the A-parameter, we may define a renormalization group invariant 
quark mass, M, by the asymptotic behavior of m, 

M = lim m(u)[2b Q g(fj,) 2 }- d °/ 2bo . (II.3.3) 

fl — >oo 

It is an easy exercise to show that M does not depend on the renormalization scheme. 
It can be computed in the SF-scheme and used afterward to obtain the running mass 
in any other scheme by inserting the proper (3- and r-functions in the renormalization 
group equations. 

To compute the scale evolution of the mass non-perturbatively, we introduce a new 
step scaling function, 

Z P (2L) 



£p(w, a/L) 



Z P (L) 



(II.3.4) 

9 2 (L)=u 



with Zp of eq. (1.1.41). Results for Sp at finite lattice spacing and the extracted 
continuum limit are displayed in Fig. 20. Details may be found in [23,98]. 
Applying crp and a recursively one then obtains the series, 

m(2- fc L max )/m(2L max ) , k = 0, 1, . . . , (II.3.5) 

up to a largest value of k, which corresponds to the smallest g that was considered in 
Fig. 20. From there on, the perturbative 2-loop approximation to the r-function and 
3-loop approximation to the /3-function (in the SF-scheme) may be used to integrate 
the renormalization group equations to infinite energy, or equivalently to g = 0. The 
result is the renormalization group invariant mass, 

~T~(g) d 



M = fn (2b g 2 )- d °/ 2bo exp l~J^9 



o l/3(g) b g 



(II.3.6) 
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Figure 20: Lattice spacing dependence of the step scaling function for the quark mass 
for Nf = 2. The coupling u ranges from u = 0.979 to u = 3.33. 



In this way, one is finally able to express the running mass fn in units of the renormal- 
ization group invariant mass, M, as shown in Fig. 21. Remember that M has the same 
value in all renormalization schemes, in contrast to the running mass fn. 

The perturbative evolution of the quark masses follows very accurately the non- 
perturbative results down to rather low energy scales. Of course, this result may not be 
generalized to running masses in other schemes. 

The point at lowest scale fj, in Fig. 21 corresponds to 

M/m = 1.296(16) at L = 2L max . (II.3.7) 

Remembering the very definition of the renormalized mass eq. (1.1.38), one can use this 
result to relate the renormalization group invariant mass and the bare current quark 
mass m on the lattice through 

M = m x 1.296(16) x Z A (g )/Z P (g , 2L max /a) = Z M (g ) m . (II.3.8) 

In this last step, one then inserts the bare current quark mass, e.g. of the strange 
quark, and extrapolates the result to the continuum limit. The bare current quark 
masses themselves are the ones for which the appropriate pseudo scalar masses are fixed 
to their experimental values. The presently available results from this strategy are listed 
in Table 4. Note that the computation of the b-quark mass required a detour through 
Heavy Quark Effective Theory, the subject of the following lecture. 

For the charm and the beauty quark masses determinations with Nf > and 
NP renormalization are still missing. However, in our opinion it is even quite early 
concerning the determinations of the light quark masses. Although some iVf -dependence 
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Figure 21: The running quark mass as a function of fx = l/L for Nf = 2. 

of the strange quark mass has been reported in the literature, one can presently not 
exclude that this is due to perturbative uncertainties or discretization errors. Note that 
the quark masses computed in the quenched approximation were in a similar stage in 
1996 but very soon afterward the uncertainties shrunk by an order of magnitude due 
to NP renormalization and continuum extrapolations. This remains to be achieved for 
the real theory with Nf > 0! 

II. 4 Renormalization scale dependence of other composite operators 

Due to our definition of the renormalized quark mass, its scale dependence is given by the 
one of the composite operator P a {x). Other composite operators can be considered and 
indeed the strategy described here has been applied to 4-fermion operators in the weak 
effective Hamiltonian [11,99-101], the HQET axial current [9,10] and in the operator 
which yields (x) of the non-singlet structure functions [125-127]. It is worth pointing 



i 


N f 


input 


Mi/GeV 


ref. 


strange 





m K ,r 


0.137(05) 


[33] 


strange 


2 


rriK,ro 


0.137(27) 


[23] 


charm 





m D ,r 


1.654(45) 


[123] 


beauty 







6.771(99) 


[124] 



Table 4: Quark masses determined with full NP renormalization and continuum limit. 
We use ro = 0.5 fm. 
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out that the non-perturbative scale dependence disagrees more (and significantly) from 
the perturbative one in some of these cases. They provide examples which emphasize 
that a fully non-perturbative renormalization is necessary to control the associated 
uncertainties. We show two examples in Fig. 22. For details we have to refer to the 
cited papers. 
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Figure 22: The running of the static-light axial current \i = 1/L for Nf = (left) [10] 
and a left-left 4-fermion operator (right) [100]. 

However, we want to emphasize again one general feature, which was noted already 
for the JVf = 2 running coupling. Significant differences between the perturbative and 
the non-perturbative result are present in cases, where the subsequent order of pertur- 
bation theory are very close to each other. The difference of 2/3-loop order running 
to the non-perturbative results can't be parameterized by the next order perturbative 
term with a reasonable coefficient. 
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Lecture III. 

Non-perturbative Heavy Quark Effective Theory 
III.l Introduction 

Similarly to the scale dependent renormalization covered in the previous lecture, the 
inclusion of heavy quarks in a lattice gauge theory simulation is a multi-scale problem. 
Such problems are always difficult and require the development of new techniques. By 
"heavy quark" we mean a quark whose mass is large compared to the intrinsic scale of 
QCD, Aqcd- If we take Aqcd ~ 500 MeV, the charm quark mass is about a factor 2 
higher, but the bottom quark mass has 



Here m,b is a quark mass defined at the scale [i = mb, but the scheme is irrelevant at 
this point; we will continue to use the symbol when this is the case. Since practical 
large volume simulations do not (yet) reach lattice spacings as small as a = 1/(5 GeV), 
we are faced with 



and a b-quark does not propagate properly on the lattice, at least when it is discretized 
with a standard relativistic QCD Lagrangian. 

For the charm quark, values am c <C 1 are achievable, but still care has to be taken. 
In Fig. 23 we show the RGI mass of the charm quark in the quenched approxima- 
tion, computed for three definitions, which differ at finite lattice spacing. While an 
extrapolation to a common continuum value is convincing, it is already clear from this 
figure, that the four times heavier bottom quark can't be treated this way. In fact, the 
main point is that the cutoff-effects become entirely non- linear in a 2 , when am is too 
large. This breakdown of the Symanzik expansion can be seen explicitly in perturba- 
tion theory [128]. In a 1-loop calculation, it has been estimated to happen around [128] 
amb ~ 1/2 or aM^ ~ 0.7. 

Various ways of coping with this problem have been proposed and investigated. Re- 
ferring the reader to reviews for other approaches [129-131], we directly turn to HQET. 
Already in 1987 Estia Eichten suggested that to describe the non-perturbative dynam- 
ics of hadrons with a single heavy quark, it is a good approximation to consider this 
quark to be static, i.e. it propagates only in time (in the rest frame of the hadron) [132]. 
This static approximation describes the correct asymptotics of bound state properties as 
mb — ► oo and corrections of order 1/rrtb can be included systematically. The expansion 
in 1/mt, is then given by an effective field theory [133-135]. It has been extended to 
transition form factors, e.g. between B- and D-mesons, assuming that also the charm 
quark can be described by the effective field theory [134, 136, 137]. This is done by 



m h rs 5 GeV ~ 10Aq C d • 




amb > 1 



(III.1.2) 
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Continuum extrapolation 




0.01 0.02 0.03 0.04 



(a/r ) 2 

Figure 23: The RGI mass of the charm quark in units of tq, evaluated with the 0(a)- 
improved formulation, as a function of the lattice spacing [123]. The physical input 
is the mass of the D s meson (and m^ro). The RGI-mass at finite lattice spacing is 
defined through the PCAC -mass, ecj. (1.2.21), of a charm-charm correlator (m cc ) or a 
charm-strange correlator (m sc ) and through the bare quark mass, m q , and eq. (1.2.24). 

considering heavy quark fields with finite velocities in the limit of large mass. We will 
here ignore this phenomenologically very interesting possibility and restrict ourselves to 
HQET at zero velocity (The formulation of HQET in Euclidean space and at non-zero 
velocity is more subtle [138]). 

In this theory considerable progress has been made recently, which we want to ex- 
plain. In particular, in order to compute 1/mb corrections, one has to know the leading 
order (including its renormalization) non-perturbatively (Sect. III. 3. 1.2). This necessi- 
tates a so-called non-perturbative matching of effective theory and QCD. A strategy for 
this [139,140] will be described in Sect. III. 5 including results for 1/mb corrections in a 
test case. The strategy solves at the same time the problem of a proper definition and 
computation of 1/mb-corrections, which is present in any regularization of the theory, 
and the problem of power divergences (~ a~ n ) in a theory with a hard cutoff such as 
the lattice regularization. After a more general introduction we will concentrate on the 
inclusion of 1/mb corrections. Because of limitations of space we mention the also very 
relevant developments in the static approximation only rather briefly. 

III. 1.1 Derivation of the classical theory 

We here go through some steps to derive the effective theory at the classical level. The 
main point is to see what assumptions have to be made and to present the explicit form 
of the Lagrangian. We follow the idea of [141], but work in Euclidean space since we 
are ultimately interested in the lattice theory. 
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In this section, we keep the dependence of the fields on the space-time coordinates 
implicit and also drop the label b on the quark field and its mass. We start from the 
Dirac-Lagrangian of a b-quark with a large mass, m, in the continuum, 

C = ^(D^ + mW (III.1.3) 
= ^Vi>, V = mj + D + j D k j k . (III.1.4) 

The light quark fields and gauge fields are not touched by our considerations. We 
write tj)\ but it is just another independent Grassmann integration variable in the path 
integral. Since we are considering the classical theory, we can assume that the fields 
are smooth. We can therefore perform an expansion in D^. More precisely, we have to 
refer to a special kinematical situation. We want to describe the dynamics of a hadron 
containing one heavy quark, where the hadron is at rest. For infinite mass, the heavy 
quark propagates only in time. Denoting the expansion parameter by e, the dynamics 
thus dictates 

D /m = 0(1) , D k /m = 0(e) , (III.1.5) 

when these derivatives act on the heavy quark fields. This is often called a power count- 
ing scheme. In the quantum theory we will have e = Aqcd/^- Obviously quantities 
such as Fuv = 0(1) are not touched by this consideration. At the lowest order in this 
expansion the ("large components") quark field (P± = ^ a ) 

rP h = P+iP, ? h = ?P+, , (m.1.6) 

propagates forward in time, while the anti-quark field, 

^ = P^, fo=^P_, (III.1.7) 

propagates backward. In a somewhat sloppy notation we will often write 0(l/m) instead 
of 0(e). The 0(l/m) terms in the Lagrangian 

C = /t at + £| tat + 0(i) (III.1.8) 

£| tat = ^ h (D + m)^ h , £f 1 = M-D + m)V E , (HI.1.9) 

connect quark and anti-quark fields. They can be decoupled through a Foldy-Wouthuysen 
rotation, 

C = $V'(j), (j) = e s ip , <fl = ^t e -S (III.1.10) 
V' = e s Ve~ s , S = ±D klk = = 0(i) , (III.l.ll) 

which yields explicitly 

V = V + ^[D klk ,V} + ^[D ni ,[D k7k ,V]] + 0(^) (III.1.12) 
= V + ^PklhM - ^[D llh70 D k7k } + O(^r) 

= 7o {70A) + m+ ^(-D k D k - ^F M a M ) + ^PfcoToTfc j + °(^) • 
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The Lagrangian then reads 

£ = £ Sh tat +£ |tat + | £ (l) +£ [l) +£ (l}| + 0( ^ ) (UL.1.13) 

4° = ^M-DkDk-^FuaM)^, (III.1.14) 
^ = \h^1u] , Fki = [D k ,Di] . (III. 1.15) 

For hadrons (or correlation functions) with a single b-quark (or anti-b-quark) only 
double insertions of C^J contribute. These are of order 1/m 2 and may be dropped at 
the order written explicitly. 

For later convenience we introduce the short hand 

4^ = -2^(°kin + Ospin), (III.1.16) 

Okin = ? h D fe -D fe ^ h =? h D 2 ^ h) (III.1.17) 

Ospin = iph-^Fkicrkiiph = ^cr-R^h- (III.1.18) 

We note that C, eq. (III. 1.13) is a low energy effective Lagrangian [142-144]. It describes 
the long wave length modes of the fields accurately and makes truncation errors, which 
are of increasing relevance for shorter wave lengths. This becomes particularly apparent 
when we remove the mass terms from the static Lagrangian and define 

£f at = ^ h (A) + e)^h , £f* = ^h(-A) + e)V>E , (III.1.19) 

where the limit e — > 0+ is to be understood in order to select the proper propagation in 
time. Replacing eq. (III. 1.9) by eq. (III. 1.19) corresponds exactly to an energy shift by 
an amount m of all states containing a single heavy quark or anti-quark. For Euclidean 
correlation functions it just leads to an additional factor of exp(— m (yo — xq)) for cor- 
relation functions where a quark propagates from xo to yo > xq. (For the anti-quark 
there is a factor exp(+m (yo — xq)) with yo < xq). 

We note again that the essential assumption is eq. (III. 1.5), namely the spatial 
covariant derivatives are counted as small compared to the mass term and the time 
derivative. This is the correct physical situation in a frame where the hadron is at rest 
and therefore at lowest order also the quark is at rest. 

Instead of carrying out the expansion of the action, one could also expand the heavy 
quark propagator in terms of 1/m. 

Quantum fluctuations are not smooth and invalidate the above "derivation" . How- 
ever, one expects that they do not modify the structure of the effective Lagrangian, but 
rather only modify the coefficients of the various terms by non-trivial renormalizations 
due to these short distance fluctuations. After all, arguing heuristically, long wavelength 
terms have been identified correctly and are described by local interaction terms. In 
local quantum field theory, also effective local quantum field theory, such terms are 
renormalized by a renormalization of the coefficients of the local fields. Below, we will 
discuss this in some detail. 
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III. 2 The effective quantum field theory 



III. 2.1 The static approximation and its symmetries 

We start with the lowest order and just consider the heavy quark; the antiquark action 
and propagator is completely analogous. In continuum Euclidean space, the classic 
Lagrangian density £| tat , eq. (III. 1.19), contains local fields of a mass dimension D < 4. 
It is power counting renormalizable. The static effective field theory is thus expected 
to be renormalizable in the usual sense, i.e. by a finite number of counter-terms. The 
possible counter-terms are restricted by the symmetries. Apart from the usual ones 
(parity, gauge symmetry . . . ), there is a well known invariance under spin rotations 
[134, 136, 137, 148]. The infinitesimal variations of these transformations can be written 
as 

<5^h = crfc^h , $ip h = -ip h a k , a k = -^e ijk aij , (HI.2.1) 
with [<Tfc, ai] = iekim&m- In addition, the action is invariant under phase transformations, 

^ _^ e ^(x) ^ ^ _^ ^ h e~^\ (III.2.2) 

with an arbitrary space- (but not time-) dependent parameter r/(x). This invariance 
corresponds to the local conservation of b-quark number, ensuring that the quark prop- 
agates only in time. The only counter-term of dimension D < 4, which involves V'hj V'h 
and is invariant under these symmetries is tph^Ph- Denoting its coefficient by 8m, the 
formal continuum quantum Lagrangian is thus simply 

£ s h tat = V^ h (A) + Smfa . (111.2.3) 

In order to discuss the equivalence of the effective theory and QCD for correlation 
functions, we introduce also the time component of the axial current in the effective 
theory, 

^ tat (x) = ^(xhoTs^x) , (ni-2.4) 

as a prototype for a composite field. There is no other operator with its quantum 
numbers and with dimension D < 3. Hence Aq \x) renormalizes multiplicatively. 



///. 2.1.1 Lattice formulation 



Just like the continuum theory we formulate the lattice theory in terms of formally 
4-component fields, satisfying 



The time-doubler is removed by choosing the lattice backward derivative 



Z) w ^ h (x) = - 



^h(x) - W\x - aO, 0)ip h (x - aO) 



(HI.2.5) 



(HI.2.6) 



61 



in the action 

S ™ = fl4 T— —x E^h(^)Po W + Srnw)Mx) ■ (m.2.7) 

1 + a 0?n w ^ 

With W^(;c,0) = U(x,0), the standard time-like links, this is the Eichten-Hill action 
[133]. For the Monte Carlo evaluation it is however of a considerable advantage to 
define the theory with more general parallel transporters W(x, 0), equivalent to U(x, 0) 
up to 0(a 2 ). In this way statistical errors of correlation functions at large separation 
xq can be reduced exponentially (in xq) and at the same time discretization errors have 
been found to be somewhat smaller [145,146]. Concentrating on conceptual issues, we 
refer to the cited papers for details on this more practical issue. 

It is an easy exercise to show that the static propagator in the presence of gauge 
fields is (aSm = ln(l + aSm)) 

Gh(x,y) = 8(x -y ) <5(x - y) exp ( - 5m w (x - y )) P w (y,x) t P + , 

(III.2.8) 

V w (x,x) = 1, V w (x,y + afi) = V w (x,y)W(y,ri, 

when the fields are normalized as in eq. (III. 2. 7). Here 0(xo),5(x) are straight forward 
lattice transcriptions of the continuum 8- and <5-functions. While in the continuum 
case, an e-prescription is necessary to select the forward propagation 0(xq — yo), the 
propagator of the lattice action with backward derivative, eq. (III. 2. 6), has this property 
automatically. The symmetries eq. (III. 2.1) and eq. (III. 2. 2) are preserved by the lattice 
regularization eq. (III. 2. 7). 

Improvement 

With these symmetries one easily goes through the steps introduced in Sect. 1.2.2 to 
find the structure of Symanzik's effective action. It turns out that the only allowed 0(a) 
terms are proportional to the light quark masses [9]. They are therefore numerically 
not very relevant and in addition their coefficients vanish for Nf = 0. It is interesting 
to note that through the static effective theory one can also give a convincing argument 
that the force between static quarks is free of linear a-effects, if the light quark action 
is 0(a)-improved [93]. 

For ease of notation we drop the sub- and superscript W from now on. 

III. 2.1.2 Renormalization 

The only term needed for the renormalization of the action, 5m, has been included above. 
The explicit form of the propagator shows that 5m enters in a purely kinematical way, 
just as an energy shift by an amount 5m compared to the unrenormalized, 5m = 
case. Thus all masses of hadrons (we consider only those with a single heavy quark) 
have the same shift and their splittings can be predicted in static approximation up to 
AQ CD /mb corrections without adjusting any parameter except for those in the "light 
part of QCD" . Obviously, the splittings are also the same if the b-quark is replaced by 
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the charm quark, but now up to Aq CD /m c , if m c is large enough for the l/m c expansion 
to make sense. 

But also in the static approximation non-trivial renormalizations occur as soon as 
hadronic matrix elements are considered. The most prominent example is the computa- 
tion of the B-meson decay constant, Fb- In QCD it can be obtained from the correlation 
function 21 

C AA (x ) = Z 2 A a 3 £ (4>(x)(4,)t(0)} (IH.2.9) 

X 

with the heavy-light axial current in QCD, = V'lT/iTs'^b) an d = VvfyTsV'i- In 
QCD, Za(<7o) is fixed by chiral Ward identities [20,21], which shows that it does not 
depend on a renormalization scale. The decay constant may e.g. be obtained from (we 
use the finite volume normalization (B\B) = 1 for the zero momentum state \B)) 

[$ QCD ] 2 = F£m B (III.2.10) 
= 2L 3 [(B\Z A A \0)] 2 = 2 lim exp(x IaaM^aa^o) , 

Xq — >QO 

where the effective mass (dof{xo) = ^[f( x o + a) — f(%o — a )]) 

r A A(xo) = -do ln(C AA (x )) (III.2.11) 

appears. In the static approximation, the QCD current is represented by 

(Ag at ) (x) = Zf*\g ^a)Ap\x) , (III.2.12) 

where we have allowed for a /U-dependence of the renormalization. This is expected 
because chiral symmetry is not just broken softly, but in the effective theory it is not 
present at all. We will come back to it in Sect. III.3.1. 

III. 2. 2 Including 1/mb corrections 

We work directly in lattice regularization. The continuum formulae are completely 
analogous. The expressions for Okim Cspim eq. (III. 1.17), eq. (III. 1.18), are discretized 
in a straight forward way, D^D^ — ► VfeVfc and F^i — > F^i with the latter defined in [19]. 
Of course other discretizations of these composite fields are possible. 

Apart from the terms in the classical Lagrangian, renormalization can in principle 
introduce new local fields compatible with the symmetries (but not eq. (III. 2.1) and 
eq. (III. 2. 2), which are broken by S pin> Okin) an d with dimension D < 5. Also the field 
equations can be used to eliminate terms. With these rules one easily finds that no new 
terms are needed and it suffices to treat the coefficients of O sp i n , £?kin as free parameters 
which depend on the bare coupling of the theory and on m^. 

21 It is technically of advantage to consider so-called smeared-smeared and local-smeared correlation 
functions, but this is irrelevant in the present discussion. 
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The 1/mb Lagrangian then reads 

4 1} (z) = -(wkinCkin(x) +U B pta0 8pta (x)) . (III.2.13) 

Since these terms are fields of dimension five, the theory defined with a path integral 
weight P oc exp(— a 4 X^tAightf^) + C^ at (x) + C^\x)]) is not renormalizable. In per- 
turbation theory, new divergences will occur at each order in the loop expansion, which 
necessitate to introduce new counter-terms. The continuum limit of the lattice theory 
will then not exist. However, that effective theory is NRQCD not HQET. Since the effec- 
tive theory is "only" supposed to reproduce the 1 /m^ expansion of the observables order 
by order in 1/mb, we expand the weight P in 1/mb, counting Wkin = 0(l/mb) = u> sp i n . 
This defines HQET. The same step has already been used in Symanzik's effective theory. 
Up to and including 0(l/mb), expectation values in HQET are defined as 



stat 

x x 
= <0) st at + Wkm(C)kin + W spin (0) spin , (III.2.14) 

where 

(O)stat = \\ O exp(-a 4 ^[Aight(x) + A? at (*)D (III.2.15) 

£ J fields x 

is defined with respect to the lowest order action, which is power counting renormaliz- 
able. The path integral defining the average extends over all fields and the normalization 
Z is fixed by (l) stat = 1- 

In order to compute matrix elements or correlation functions in the effective theory, 
we also need the effective composite fields. At the classical level they can again be 
obtained from the Foldy-Wouthuysen rotation. In the quantum theory one adds all 
local fields with the proper quantum numbers and dimensions. For example the effective 
axial current (time component) is given by 

A* QET (x) = Z^ QET [^ tat (x) + C A QET <5Ao at (^)], (III.2.16) 
SA^(x) = ?i(x)^(V l + V*) 7l 75^h(x). (III.2.17) 

Before entering into more details on the renormalization, we show some examples how 
the 1/mb-expansion works. 

III. 2. 2.1 1/mb- expansion of correlation functions and matrix elements 
For now we assume that the coefficients 

0(1): 6m, Z^, 

(III.2.18) 

0(l/m b ) : ^hqet 
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are known as a function of the bare coupling g$ and the quark mass m^. Their non- 
perturbative determination will be discussed later. 

The rules of the 1/mb-expansion are illustrated for of Caa(xo), eq. (III. 2. 9). One 
uses eq. (III. 2. 14) and the HQET representation of the composite field eq. (III. 2. 16). 
Then the expectation value is expanded consistently in 1/rat,, counting powers of l/m^ 
as in eq. (III. 2. 18). At order l/m^, terms proportional to oj k i n x c a ^ ET etc. are to be 
dropped. As a last step, we have to take the energy shift between HQET and QCD 
into account. Therefore the correlation function obtains an extra factor exp(— xom^), 
where the scheme dependence of is compensated by a corresponding one in 8m. One 
arrives at the expansion 

Caa(xo) = e-~(Z™) 2 [ci t l t (x ) + ^ QET CI A i A (x ) (111.2.19) 

+ UJ kin C A a(x ) + ^spin CaT^o) 

with (remember the definitions in eq. (III. 2. 14)) 

cffisxo) = << at ( a 0(*< at (o)) + W + (^(aOC^WW, 

Cfl{x Q ) = (^ tat (x)(^ tat (0))t) kin , (III.2.20) 



C A t(x ) = {At*\x){At*\W) 



spin • 

22 



It is now a straight forward exercise to obtain the expansion of the B-meson mass 

m B = m h + 5m + E stat + Wkin^kin + Wspin^spin , (III. 2. 21) 



E stat = - lim d lnCfl\x ) 

XQ — >OC 



(III.2.22) 



6m=0 

E ^ = " ii^^Pkin^o), Pkin(xo) = *&) 0) , (IIL2.23) 
£ spin = -^aopspin^o), Pspin(xo) = , (III.2.24) 



and its decay constant 



FbV^b = lim {2exp(m B x )C AA (x )} 1/2 , (III.2.25) 



^1/2 

= zj* QET Jm^ $ stat { 1 + \x G [uj kin E kin + uj spin E kin } 

+ 2 C A P<5a(^o) + 2 CJ kinPkin(a;o) + 2 ^spinPspin^o)} , (III. 2. 26) 

$ " at = l™ oo {2exp(ii; stat ,o)CIt(xo)} 1/2 , P6A(x ) = ^§^. 



22 It follows from the simple form of the static propagator that there is no dependence on 8m except 
for the explicitly shown energy shift 8m. 
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Using the transfer matrix formalism (with normalization {B\B) = 1), one further ob- 
serves that 

E khl = -{B\o? CWO, z)|B)rtat, 

(III.2.27) 

E spin = -(B\a 3 J2°s P in(0,z)\B) stat . 

z 

As expected, only the parameters of the action are relevant in the expansion of hadron 
masses. 

A correct split of the terms in eq. (III. 2. 21) and eq. (III. 2. 26) into leading order 
and next to leading order pieces which separately have a continuum limit requires more 
thought on the renormalization of the 1/mb-expansion. We turn to this now. 

111.2.2.2 Renormalization and continuum limit 

For illustration we first check the self consistency of eq. (III. 2. 19). The relevant question 
concerns renormalization, namely: are the "free" parameters 5m . . . c^ ET sufficient to 
absorb all divergences on the r.h.s.? We consider the most difficult term, C^™(xo). 
According to the standard rules, it is renormalized as 

fe) R (* ) = (^I tat ) 2 x (III.2.28) 



a 7 £ (a^(x) «* (0))t (Okta) R W) + C.T 



where C.T. denotes contact terms to be discussed shortly. The operator (O k i n ) R (z) 
involves a subtraction of lower dimensional ones, 

(O kin ) R (z) = Z 0kin (O kin (z) + ^ h (z)D ip h (z) + ^J\(z)Mz)) , 

(III.2.29) 

written here in terms of dimensionless q. Since we are interested in on-shell observables 
(xq > in eq. (III. 2. 19)), we may use the equation of motion Dqi/j^z) = to eliminate 
the second term. The third one, ^ijj- ll (z)ijjh(z), is equivalent to a mass shift and only 
changes 5m, which is hence quadratically divergent 23 . Thus all terms which are needed 
for the renormalization of Okin are present in eq. (III. 2. 19). 

It remains to consider the contact terms in eq. (III. 2. 28). They originate from sin- 
gularities in the operator products O k x n {z)A S Q a,t {x) as z ~~ * x ( an ^ Okm{ z )(A ! Q a ' t )^ {0) as 
z — > 0), in complete analogy to the discussion in Sect. 1.2.2.3. Using the operator prod- 
uct expansion they can be represented as linear combinations of A l Q a ' t (x) and SAq 1 (x). 



23 Using the explicit form of the static propagator, eq. (III. 2. 8), one can check that indeed 
« 3 Ex (^o tat (^)(^o tat (0)) t a 4 E z V ; h( 2 )'/'h( 2 )) = XoCJ^fxo), which can be absorbed by a l/m b 

\ / stat 

correction to 5m. 
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Such terms are contained in eq. (III. 2. 19) in the form of an d C|aa- Indeed Ag tat (x) 
and <5^4Q tat (x) are the only operators of dimension 3 and 4 with the correct quantum 
numbers. Higher dimensional operators contribute only terms of order a. 

Note that the coefficient of A l Q a,t (x) in the expansion of the operator product 
a 4L ^2 z O\an{z)A S Q &t {x) is power divergent ~ 1/a, for simple dimensional reasons. This 
means that there is a power divergent contribution to Z A . As this happens only at 
order l/m^, not at the lowest order, this contribution to Z A behaves like ~ l/(arab) 
for small lattice spacing. 

We conclude that all terms which are needed for the renormalization of C A ™(xo) 
are present in eq. (III. 2. 19); the parameters may thus be adjusted to absorb all infinities 
and with properly chosen coefficients the continuum limit of the r.h.s. is expected to 
exist. The basic assumption of the effective field theory is that once the finite parts of 
the coefficients have been determined by matching a set of observables to QCD, these 
coefficients are applicable to any other observables. 



111.2.2.3 The flavor currents in the effective theory 

For later use we here give the expressions for the heavy-light currents. They are relevant 
in weak B-meson decays. For better readability, we include again the time component 
of the axial current. 

Following our general rules for finding the HQET fields which represent the QCD 
ones we find 

Aq QET (x) = ZA QET [A s tat (x)+CA QET ^ tat (x)], (III.2.30) 
y HQET (x) = Z^ QET [T/ stat (x) + C y QET 5T/ stat (x)], (III.2.31) 
vf QET (x) = Z^ QET [T/ fe stat (x) + 4 QET <5T4 stat (x)], (III.2.32) 
Aj QET (x) = Z^ QET [A| tat (x) + 4 QET M| tat ( 2; )]. (III.2.33) 

Ignoring that one can in principle simplify, for instance ip\ (a:) 7075^ ( x ) = ~' t Pi( x )l5' l l J h( x ) 
due to P+iph = V'hj the basis fields are written in full analogy to the ones in QCD: 

V£ ut = ?i7oV<h, ^ s tat = ?i7075^h , (III.2.34) 
Ff at = ? l7 ^h, 4^=^7*75^, (III.2.35) 
5V ^ = ^^7* 7 .^ h) M stat = ^^ i i^* 7 . 75V;h) (in.2.36) 

<5^ stat = -^^±^ 7i7fc ^, M| tat = V; i %^ 7 7 fc 75^h, (IH.2.37) 
We have chosen the bare fields such that they are exactly related by spin rotations, 24 

rfcrT/stati • a stat x^rxT/stati -r /istat 



24 The unnatural — sign in eq. (III. 2. 37) is present because we remain with the definition of SAq^ 
in [9] and do not want to introduce signs in eq. (III. 2. 38). 
The other rotations look like SiAf^ = i e fcj i V, stat - i SkjAf^. 
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(III.2.38) 

^fc^statj = -^stat^ S*[6A$ at ] =i(5V fc Stat . 

In general, all Z-factors and coefficients c H( ^ ET are functions of go and am^ which are 
to be determined by matching non-perturbatively to QCD. Eq. (III. 2. 38) will be needed 
in Sect. III. 3.1 when we discuss the static limit of the currents. 

III. 2. 3 Schrodinger functional correlation functions 

For an understanding of the details of the tests of HQET (Sect. III. 4) as well as the 
non-perturbative matching to QCD (Sect. III. 5) we will also need some Schrodinger 
functional correlation functions and their HQET expansion, which have not been defined 
yet. We give these details [140] now. The reader who is only interested in the general 
concepts may skip this section. 

In [9] static quarks in the Schrodinger functional were discussed including Symanzik 
0(a)-improvement. It turns out that there are no dimension four composite fields which 
involve static quarks fields and which are compatible with the symmetries of the static 
action and the Schrodinger functional boundary conditions and which do not vanish by 
the equations of motion. Thus there are no 0(a) boundary counter terms with static 
quark fields. For the same reason there are also no 0(l/mb) boundary terms in HQET. 
This then means the HQET expansion of the boundary quark fields £ is trivial up to 
and including 1/rrib terms. 

In the spatial boundary condition of the fermion fields, 

ip(x + kL) = e i6 tp(x) , ijj(x + kL) = e -<fl ^(aj) , (III.2.39) 

the same phase is taken for all quark fields, whether relativistic or described by 
HQET. 25 In QCD, relevant correlation functions in the pseudo-scalar and vector channel 
are 

/a(*o,0) = -yE ((^i)o(*Kb(y)75Ci(z)) , (III.2.40) 
k Y (x ,9) = ~ J2 ((W*Kb(y)7feCi(z)) , (III.2.41) 

with the 0(a) improved currents [139] 

(Ai)fj, = tpa^ipb + acAd^tpi^iph (III. 2.42) 

{ydn = V'lTAiV'b + acyd^iia^ipi, . (III. 2. 43) 

25 In principle we can easily have different 9 for different fields, as long as they are quenched, so in 
particular for the heavy field. Since this freedom has not yet been used, we do not discuss it here. In 
applications with dynamical fermions it may well become relevant. 
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Furthermore we consider boundary to boundary correlation functions 

2L 6 



/iW = E (Ci / (u) 7 5CL(v)C b (y)75Ci(z)) , (IH.2.44) 



u,v,y,z 

h(9) = ~ Yl (Ci'(u) 7fc CUv)C b (y)7 fc Ci(z)) . (III.2.45) 

u,v,y,z,fc 

Their renormalization is standard [46], for example at vanishing light quark masses 

[fA] R (x ,9) = Z A (l + i6 A am qib )Z c 2 (l + 6 c am qib )/ A (xo^), (III.2.46) 
] R (6) = Z ? 4 (l + 6 c am qj , 



[A] R (0) = ZHl + 6 c am q , b ) 2 /iW, (III.2.47) 



with Zq a renormalization factor of the relativistic boundary quark fields and another 
improvement coefficient. 

Their expansions to first order in l/m b read 

[/ A ] R = Zf^Z^e-^^f^ + ^f^ + u^f^ + LJ^f^} 

(III.2.48) 

[Mr = Z™Z ?h Z ? e-™ {a^ + 4 QET k^ + u; km k^ + u spin kf n } 

(III.2.49) 

V HQET ry ry — m h xn r pstat . HQET ^stat . ^kin 

= -Z y ^ Z (h Z c e b °{/ A +c v ^ /,5 A +w kin / A 

— h i ■ f s P in \ 

[/i] R = Z 2 h Z 2 e- m ^{/f at +a; ki n/i kin + ^pin/r in } , (III.2.50) 

[*i] R = ^Zle-^^j/f^+^/^-i^pin/r 11 } • (III.2.51) 



Apart from 



/ 4 r(x o ,0) = 4E (^o at (^)Ch(y)75Ci(z)) (III.2.52) 



2 



the labeling of the different terms follows directly the one introduced in eq. (III. 2. 14). 
We have used identities such as /^ in = — ky n , / A pm = 3fcy in . As a simple consequence 
of the spin symmetry of the static action, these are valid at any lattice spacing. 



III. 3 The scope of the theory 

We are now in the position to discuss what can be done in the effective theory. This 
concerns also the continuum effective theory and in particular the question, where per- 
turbation theory is sufficient to give an answer and where it is not. 
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III. 3.1 A first example: the decay constant 



///. 3.1.1 Renormalization and matching in perturbation theory 

The matrix element $Q CD , eq. (III. 2. 10), is scale independent, due to the chiral symme- 
try of QCD in the massless limit. But of course it depends on the mass of the b-quark. 
In the effective theory this symmetry is absent and, remaining a while at the lowest 
order in l/m b , Z^ at depends on the renormalization scale, /U, used in the renormal- 
ization condition which defines the finite current. On the other hand, as long as the 
effective theory is considered by itself, Z^ tat does not depend on the mass of the b-quark 
in an obvious way. The mass-dependence comes in by choosing an appropriate finite 
renormalization such that the matrix elements of the current in the effective theory are 
equal to the ones of the QCD current (up to 0(l/mb)). This step is called matching. 
Choosing to renormalize the current in the effective theory in the MS scheme (any other 
scheme would of course be possible) we therefore have 



$Q CD = C match (m b ,/i) x + 0(l/m b ) 



with 



a 



match 



(m b ,/i) = 1 + ci(m b //i)g|g(/i) + . . . 



(III.3.1) 



(III.3.2) 



The finite renormalization factor C ma t c h is determined (usually in perturbation theory) 
such that eq. (III. 3.1) holds for some particular matrix element of the current and will 
then be valid for all matrix elements. Obviously, the //-dependence in eq. (III. 3.1) is 
artificial, since we have a scale-independent quantity in QCD. Only the m b -dependence 
is for real. 

The /i-dependence in eq. (III. 3.1) is removed explicitly by changing from $j^g(/i) = 



2 , ^ ta lg(0|^o|-B) stat to the RGI matrix element 



<I> 



RGI 



lim 

fi— »oo 



2b g 2 (fi) 



-70/260 



^Ms(^) • 



Here, the lowest order coefficient of the /3-function, bo enters as well as 70 
defined by 



lid) 



vstat 



We can now write down the HQET-expansion of the QCD matrix element 

$Qcd = Cps (M h /A m ) x <D RG i + 0(l/m b ) . 
The relation between eq. (III. 3. 5) and eq. (III. 3.1) is easily seen by using 



$ MS(^) 



2&oS 2 (m) 



-7o/2feo 



cxp 



7Ms(g) _ 7^ 



(III.3.3) 

-l/(4vr 2 ) 

(III.3.4) 
(III.3.5) 
(III.3.6) 
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setting the arbitrary renormalization point fi to and identifying 



Cps(^) = C match (D%^ (IH.3.7) 
\ a ms/ *RGi 



2b g 2 (m b ) 



7o/2&o 

exp 




7match(g) _ _70_ 

^ms(5) b g 



where g is taken in the MS scheme. The last equation provides a definition of the anoma- 
lous dimension 7 ma t c h m the "matching scheme". Perturbatively, it has contributions 
from Tj^g as well as from C ma t c h> namely 

7match(s) = -JOS 2 ~ [7P + 2W1)]5 4 + • • • • ( IIL3 - 8 ) 

Replacing the MS coupling by a non-perturbative one, 7match may also be defined beyond 
perturbation theory through eqs. (III. 3. 7, III. 3. 5). 26 Another advantage of eq. (III. 3. 7) 
(compared to eq. (III.3.1)) is that Cps is independent of the arbitrary choice of renor- 
malization scheme for the composite operators in the effective theory. Apart from the 
choice of the QCD coupling, the "convergence" of the series eq. (III. 3. 8) is dictated by 
the physics, nothing else. 

Note further that (at leading order in l/mt>) the conversion function Cps contains 
the full (logarithmic) mass-dependence. The non-perturbative effective theory matrix 
elements, ^rgp axe mass independent numbers. Conversion functions such as Cps are 
universal for all (low energy) matrix elements of their associated operator. Thus 



Caa(xo) X0>> ~ /mh [Cps(^)n t , a RGi] 2 ^o tat (^) t ^o tat (0)) + 0(^) 

MS (III.3. 



is a straight forward generalization of eq. (III. 3. 5). Here, 

Zf% Gl = lim^oo [26 5 2 (m)]" 70/26 ° 2*gg(/*). 

Analogous expressions for the conversion functions are valid for the time component 
of the axial current replaced by other composite fields, for example the space compo- 
nents of the vector current. Based on the work of [147-149] and recent efforts their 
perturbative expansion is known including the 3-loop anomalous dimension 7match 

ob- 
tained from the 3-loop anomalous dimension 7^ [150] and the 2-loop matching function 
Cmatch [151-153]. Figure 24, taken from [154], illustrates that the remaining 0(^ 6 (mb)) 
errors in Cps seem to be relatively small. 

We return to the full set of heavy-light flavor currents Sect. III. 2. 2.3. The bare 
fields satisfy the symmetry relations eq. (III. 2. 38). The same is then true for the RGI 
fields in static approximation. Furthermore, the axial currents are related to the vector 



Clearly the r.h.s. of eq. (III. 3. 7) is a function of g (m-b), i.e. a function of mb/Aj^g. We prefer to 
write it as a function of the ratio of renormalization group invariants, Mb/A^g. 
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Figure 24: Cps estimated in perturbation theory. For B-physics we have Aj-^/Mb ~ 0.04. 



ones by a chiral rotation of the light quark fields [155]. It then follows that in static 
approximation the effective currents are given by 

^hqet = Cps (M h /A m )Zf% G1 (g )A^, (III.3.10) 

V^ ET = C v (M h /A Ws )Zt% Gl (go)Vt a \ (IH.3.11) 

^HQET = C y(M h /A m ) Zf$ Gl (g ) Z$? A (g ) V*« , (III.3.12) 

A* QET = C A (M h /A m )Zf$ G1 (g )Z^ A (g )ArK (III.3.13) 

The finite renormalization ZyJ^(go) can either be computed by a Ward identity [155] or, 
if one has a regularization with exact chiral symmetry, it is one. The factor Z\ tl R G i(go) 
is analogous to Zm(9o) in eq. (II. 3. 8) and has been computed in a similar way [10]. It 
can be split into 



^Gl(2o) = |gy x Zf*\g ,2L max /a) 



(III.3.14) 

At=l/(2L max ) 



where the second factor depends on the lattice action, the first one does not. The first 
factor has been shown in Fig. 22 for a series of \x in the quenched approximation with 
a smallest scale of fj, = l/(2L max ). 

We finally note that the bare pseudo scalar density, P, and the scalar density, S, 
are identical to — ^4o tat an d Vo tat , respectively. This then also holds for the RGI ones; 
they are given by -Z^ GI ( 5o )^ s tat and Zf% Gl {g ) Z^f^V^ . 

For the renormalized _p H Q ET and 5 H< ^ ET , there is in principle an arbitrariness, since 
in QCD they depend on a scale and a scheme. This arbitrariness is fixed by considering 
the RGI fields in QCD. They satisfy the PCAC and PCVC relation with the RGI quark 
masses, for example (Mi is the mass of the light quark) 

d^An)^ = (M b + M!)P RGI . (III.3.15) 

Taking the vacuum-to- B-meson matrix element of this relation (in static approximation) 
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Figure 25: Example of an interpolation between a static result and results with rah < rat,. Continuum 
extrapolations are done before the interpolation [156]. The point at 1/roraps = is given by t"q^ 2 $rgi- 

fixes 

pHQET = _ Cps(Mb/A _ ) ^ Z stat Gi(5oM stat ) (m 3 lg) 

5 hqet = C y(M h /A m )^Z$f A (g )Zf$ G1 (g )V^ (III.3.17) 

in static approximation. (Choosing a different matrix element would change m-Q — > 
ttib + O(Aqcd)-) We leave it as an exercise to verify these relations. 

All conversion functions C are known up to (presumably small) a(mt,) 3 errors. At 
least when the conversion functions are known with such high precision, the knowledge 
of the leading term in expansions such as eq. (III. 3. 5) is very useful to constrain the 
large mass behavior of QCD observables, computed on the lattice with unphysical quark 
masses mj, < nib, typically nih as m c h ar m- As illustrated in Fig. 25, one can then, with 
a reasonable smoothness assumption, interpolate to the physical point. 

The relation between the RGI fields and the bare fields has so far been obtained 

for 

• z a%gi with both N f = 1 10 ] and N f = 2 1 157 ]; 

• the parity violating AB = 2 four fermion operators [11, 158] for Nf = while 
Nf = 2 has been started. 

In the second case, the matrix elements of the two operators, evaluated in twisted mass 
QCD will give the standard model B-parameter for B-B mixing. 

Thus, soon one will be able to do interpolations such as Fig. 25 also for semi-leptonic 
decays such as B — ► tt I v (close to zero recoil) and for the B-parameter. 

III. 3. 1.2 Beyond the leading order: the need for non-perturbative conversion functions 

Still, getting the continuum extrapolations of the data at finite heavy quark masses in 
Fig. 25 under control may represent a challenge when dynamical fermions are included. 
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Furthermore one should not forget that the functional form of the interpolation in that 
figure does essentially assume that the 1 jm expansion remains reasonably accurate also 
significantly below m^. 

It is therefore natural to try to compute the 1/mb correction directly in HQET. 
However, if one wants to do this consistently, the leading order conversion functions such 
as Cps have to be known non-perturbatively. This general problem in the determination 
of power corrections in QCD is seen in the following way. Consider the error made in 
eq. (III. 3. 7) (or eq. (III. 3.1)) when the anomalous dimension has been computed at I 
loops and C ma t c h at I — 1 loop order. The conversion function Cps is then known up to 
an error 



As mb is made large, this perturbative error becomes dominant over the power correction 
one wants to determine. Taking a perturbative conversion function and adding power 
corrections to the leading order effective theory is thus a phenomenological approach, 
where one assumes that the coefficient of the {g 2 (rnb)] 1 term (as well as higher order 
ones) is small, such that the A/mb corrections dominate over a certain mass interval. 
In such a phenomenological determination of a power correction, its size depends on 
the order of perturbation theory considered. A theoretically consistent evaluation of 
power corrections requires a fully non-perturbative formulation of the theory including 
a non-perturbative matching to QCD, see Sect. III. 5. 

Note for experts. Eq. (III. 3. 18) has little to do with renormalons. Rather it is due 
to the truncation of perturbation theory at a fixed order. Of course a renormalon-like 
growth of the perturbative coefficients does not help. 

III. 3. 1.3 Splitting leading order (LO) and next to leading order (NLO) 

We just learned that the very definition of a NLO correction to means to take 
eq. (III. 2. 26) with all coefficients Z E( ^ ET . . . c^ ET determined non-perturbatively. We 
want to briefly explain that the split between LO and NLO is not unique. This is fully 
analogous to the case of standard perturbation theory in a, where the split depends 
on the renormalization scheme used, or better on the experimental observable used to 
determine a in the first place. 

Consider the lowest order. The only coefficient needed in eq. (III. 2. 26) is then 
Cps^a rgi = nm m b ^oo Z E< *^ ET . It has to be fixed by matching some matrix element of 
^stat ^ ma t r j x element of A$ in QCD. For example one may choose (B'\Al\0), with 
\B') denoting the first pseudo-scalar excited state. Or one may take a finite volume 
matrix element (B(L)\Aq\Q(L)), see Sect. III. 4. Since the matching involves the QCD 
matrix element, there are higher order in 1/mb "pieces" in these equations. There is no 
reason for them to be independent of the particular matrix element. So from matching 
condition to matching condition, Cps^arqi differs by ©(AQCD/^b) terms. 




(III.3.18) 
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The matrix element in static approximation inherits this 0(AQco/^b) ambigu- 
ity. The A.QCT>/ m b corrections are hence not unique. Fixing a matching condition, the 
leading order Fb as well as the one including the corrections can be computed and has a 
continuum limit. Their difference can be defined as the 1/rrtb correction. However, what 
matters is not the ambiguous NLO term, but the fact that the uncertainty is reduced 
from 0(AQCD/^b) m the LO term to 0(AQ CD /m^) in the sum. 

III. 3. 2 A second example: mass formulae 

Let us represent the formula for the meson mass, eq. (III. 2. 21), including also the vector 
meson as 

= \i m B + 3m B »] = m b + Sm + E stat + ^ kin £ kin + 0(l/mg) (III.3.19) 

and 

Am B = m B * - m B = uj spin E spin + 0(1 /ml) , (III. 3. 20) 

where the fact that C S pin does not contribute in eq. (III. 3. 19) is a consequence of the 
exact spin symmetry of the lowest order. It is a simple exercise to verify the above 
equations explicitly by considering correlation functions of Vk = V'lTfcV'b instead of Aq. 
Usually these equations are written as [159] 

_ 1 2 

m| v ~ m b + A + Ai , Am B A 2 . (III.3.21) 

2mb w-b 

We arrive at a similar form by defining power divergent subtractions for E\^ n and E stajt 
vs. 

m | v = m h + m ^ + m p + m g a) + m g 6) + 0(A| CD /m£) , (III.3.22) 



where we have split up as 



O(Aqcd) : m^=Sm 8Ut + E^ t , m|f = E stat - E^ t , (III.3.23) 

0(A^ CD /m b ) : m£ a) =6m (l) +iu kin E^, m£ b) = co kin [E kin - E^] . 

(III.3.24) 

The subtraction terms (the Hamiltonian should be defined as H = — ^ln(T) in terms 
of the transfer matrix T) 

= (/3|H|/3) stat | 5m=0 , E^ = tf\ - a 3 £ O kin (x)|/?) stat . (III.3.25) 

X 

are chosen such that ^ . . . m^ are finite and have a continuum limit. This follows 

-r^-stat 

from our discussion of renormalization. By 5m we denote 5m in static approximation 
and by 5m ^ the piece that is to be added when w k i n / 0; it accounts for the term c-ij 'a 2 
in eq. (III.2.29). 
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This rewriting exposes that the split up of the meson formulae into various orders 
of 1/mb is not unique. In the chosen form it depends on the arbitrary state \(3). In 
other words: non-perturbatively, A and Ai can be defined exactly, but in many ways. 27 

Toward the end of this lecture we will see how the subtraction terms and m b + 8m 
can be defined such that they are computable in practice and the remaining error is 
reduced to 0(AQ CD /m b ). 

III.4 Non-perturbative tests of HQET 

Although it is generally accepted that HQET is an effective theory of QCD, tests of this 
equivalence are rare and mostly based on phenomenological analysis of experimental 
results. A pure theory test can be performed if QCD including a heavy enough quark 
can be simulated on the lattice at lattice spacings which are small enough to be able to 
take the continuum limit. This has recently been achieved [154] and will be summarized 
below. We start with the QCD side of such a test. Lattice spacings such that am b <S 1 
can be reached if one puts the theory in a finite volume, I? x T with L, T not too large. 
We shall use T = L. For various practical reasons, Schrodinger functional boundary 
conditions are chosen. Equivalent boundary conditions are imposed in the effective 
theory. We then consider correlation functions such as /a and f\. The first one is the 
correlator of boundary quark fields £ (located at xq = 0) and the time component of 
the axial current in the bulk (0 < xq < T). The second one describes the propagation 
of a quark-antiquark pair from the xq = boundary to the xq = T boundary. See 
Sect. III.2.3 for details. 

We then take a ratio for which the renormalization factors of the boundary fields 
cancel, 



y PS (L,M b ) =. z a /a(l/2) 



1: 



ML)\A \B(L)) 



T=L 



MD)\\ \\\B(L)) 



\B(L)) = e-W\<pB(L)) , \n(L)) = e-W\<po(L)) . 

As shown in the above equations, Yps can be represented as a matrix element of the axial 
current between a normalized state \B(L)) with the quantum numbers of a B-meson 
and \£l(L)) which has vacuum quantum numbers. The time evolution e -iH / 2 ensures 
that both of these states are dominated by energy eigenstates with energies around 2/L 
and less (above m b ). In other words, HQET is applicable if 1/L <C m b (and of course 
Aqcd < rn h ). 

One then expects (for fixed LAqcd) 

F PS (L, M b )/C PS (M b /A) = X RG i + 0(1/ z) , z = M h L , (III.4.2) 



27 In phenomenological applications [160], one works in dimensional regularization where the subtrac- 
tion terms can be omitted. One defines a perturbative scheme for mt>. It is then clear that (for example) 

2 £-1-1 m b — 

A depends on the order of perturbation theory as [g (mb)] mb] S> Aqcd as in eq. (III. 3. 18). Quan- 
tities such as A are then effective parameters, not directly related to the asymptotic l/m,b expansion. 
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Figure 26: Testing eq. (III. 4. 2) through numerical simulations in the quenched approximation and for 
L ~ 0.2 fm [154]. The physical mass of the b-quark corresponds to z ~ 5. 



where the 1 /rrih corrections are written in the dimensionless variable 1/z and Argi is de- 
fined as Yps but at lowest order in the effective theory and normalized as in eq. (III. 3. 3). 
Of course such relations are expected after the continuum limit of both sides has been 
taken separately. For the case of Yps(L, Mb), this is done by the following steps: 

• Fix a value uq for the renormalized coupling g 2 (L) (in the Schrodinger functional 
scheme) at vanishing quark mass. In [154] uq is chosen such that L ~ 0.2 fm. 

• For a given resolution L/a, determine the bare coupling from the condition g 2 (L) = 
uq. This can easily be done since the relation between bare and renormalized 
coupling is known [98]. 

• Fix the bare quark mass m q of the heavy quark such that LM = z using the 
known renormalization factors Zm, Z in M = ZmZ (1 + a6 m m q ) m q , where Z was 
introduced in eq. (1.2.24) and Zm in eq. (II. 3. 8). 

• Evaluate Yps and repeat for better resolution a/L. 

• Extrapolate to the continuum as shown in Fig. 26, left. 

As can be seen in the figure, the continuum extrapolation becomes more difficult as the 
mass of the heavy quark is increased and 0((ara) 2 ) discretization errors become more 
and more important. In contrast the continuum extrapolation in the static effective 
theory (Fig. 27) is much easier (once the renormalization factor relating bare current 
and RGI current is known [10]). After the continuum limit has been taken, the finite 
mass QCD observable Yp$(L, M) turns smoothly into the prediction from the effective 
theory as illustrated in the r.h.s. of Fig. 26. Several such successful tests were performed 
in [154], two of them with the static (M — ► 00) limit known from the spin symmetry 
of HQET. For lack of space we do not show more examples but only note that the 
coefficient of the l/z n terms in fits to the finite mass results together with the static 
limit are roughly of order unity. 
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Figure 27: Continuum extrapolation of Xrxji [154]. 



Of course, finite mass lattice QCD results have been compared to the static limit 
over the years, see for example [161-173] and references therein. So what is new in 
the tests just discussed? The composite fields were renormalized non-perturbatively 
throughout and, by considering a small volume, the continuum limit could be taken at 
large quark masses. 



III. 5 Strategy for non-perturbative matching 

Following Sect. III. 3. 1.2, the missing piece for a general computation including 1/mb 
corrections is a practical strategy for determining the parameters in the Lagrangian 
and in the effective fields beyond perturbation theory. Let us denote the number of 
parameters which have to be determined, not counting the parameters in the light 
sector of QCD, by iV e ff. For instance, including 1/mb terms but not considering matrix 
elements of any composite fields, we have N e g = 3, namely Wkim^spin and + 5m. 
Given these three parameters, all masses can be computed. If in addition we want to 
compute matrix elements of Aq, such as Fb, we have iV en - = 5, since also Z^ ET , c^ ET 
are parameters of the theory. 



III. 5.1 Matching in small volume 

Observables, e.g. dimensionless renormalized correlation functions or energies are de- 
noted by $ H Q ET in the effective theory and by $Q CD in QCD. The N e g unknown param- 
eters can be determined from $u = 3>i7 , k = 1 . . . iV e g , provided the sensitivity 
of these conditions to the desired parameters is sufficient. In general the determination 
of ( I ) ^ CD will be very difficult because a b-quark has to be simulated and 0((ami,) 2 ) 
cutoff effects will be large. The way around is once again to consider a finite volume, 
where small lattice spacings are accessible. In practice using furthermore Schrodinger 
functional boundary conditions is a good idea, since then the simulations can easily 
be done also with dynamical fermions. Given the experience of the tests of HQET 
(Sect. III.4), a good choice is L = L x « 0.4 fm. Then 1/z = l/(LM h ) « 1/10 and the 
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0.4 fm 




| | | | | | 1 1 1 1 1 1 1 1 1 1 -I— I— | matching 

iNNNNN; f\ fftttt 

$* QET (L 2 ,M b ) $HQET (LiiMb) 

Figure 28: The strategy for a non-perturbative determination of the HQET-parameters from QCD 
simulations in a small volume. Steps indicated by arrows are to be repeated at smaller lattice spacings 
to reach a continuum limit. 

HQET expansion is very accurate. So for the matching step we impose 

^ QET (L 1 ,M h ) = ^ CD (Li,M h ), k = l,...,N eS . (III.5.1) 

to determine the iV e ff parameters in the effective theory (right hand side of Fig. 28). We 
assume that the observables &k(L, Mb) have been made dimensionless by multiplication 
with appropriate powers of L. They should be chosen with care (e.g. no large momenta 
should appear) but the effect of variations in the matching conditions on the final results 
is in any case of a higher order in the 1/mb expansion. 



HQET 

m B 
F B 

(large 
volume) 



III. 5. 2 Step scaling functions 

The matching conditions, eq. (III. 5.1), define the HQET parameters for any value of 
the lattice spacing (or equivalently bare coupling). In practice, for L\ ~ 0.4 fm, the 
parameters of the effective theory are then determined at rather small lattice spacings in 
a range of a ~ 0.02 fm to a ~ 0.05 fm. Large volumes as they are needed to compute the 
physical mass spectrum or matrix elements then require very large lattices (L/a > 50). 
A further step is needed to bridge the gap to practicable lattice spacings. A well-defined 
procedure is as follows (bottom part of Fig. 28). We define step scaling functions [37], 
F k , by 

<&* QET ( S L, M) = F k ({$J QET (L, M) , j = 1 . . . N cS }) , k = 1 . . . iV eff , (III.5.2) 
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where usually one uses scale changes of s = 2. These dimensionless functions describe 
the change of the complete set of observables } under a scaling of L — ► sL. In 

order to compute them one 

i selects a lattice with a certain resolution a/L. 

ii The specification of <I> • (L, M), j = 1, . . . , iV n , then fixes all (bare) parameters 
of the theory. 

iii The l.h.s. of eq. (III. 5. 2) is now computed, keeping the bare parameters fixed while 
changing L/a — ► V /a = sL/a. 

iv The values for the continuum Fk are reached by extrapolating the resulting lattice 
numbers to a/L — > 0. 

Starting from L = L\ it turns out that a single step going to L2 = 2L\ is sufficient [140]. 
Then one can switch at fixed bare parameters to a large volume where finite size effects 
are negligible (left part of Fig. 28). This is done in full analogy to the steps i. . . iv above. 
One only computes the large volume quantities for the bare parameters fixed in step ii 
(with L = L<i). Again a continuum extrapolation can be carried out. 

III. 5. 3 Example: The mass of the b quark 

For illustration purposes we consider a simple example, the computation of the b-quark 
mass, starting from the observed (spin averaged) B-meson mass. For this calculation one 
obviously has to consider a range of masses Mb in eq. (III. 5.1) and determine the physical 
point from the requirement mg = m^\ experiment- It is thus the first computation to be 
carried out. Subsequently, one may directly choose the physical point in eq. (III. 5.1). 

III. 5. 3.1 Static approximation 

Remembering eq. (III. 3. 23), we are after a precise definition and calculation of E^ t 
and mg = 5m + E^t- Already in the static approximation, a non-perturbative 
matching is required, if one wants to take the continuum limit. (Perturbatively one de- 
termines 5m with an in the continuum limit divergent error term of order A(5m ) = 



In eq. (III. 5.1) we have the simple case N e g = 1. We omit the discussion of fixing 
the bare light quark masses and coupling. Obviously any finite volume energy in the 
b-sector, denoted by T, will do to fix 5m. Two precise definitions are given in the 
following four equations. The reader who is only interested in the general concept may 
skip this detail. 

The first choice which comes to mind is [60, 139, 140] 



ci+m 




F av (L,0 o ) 



do + d* 



-FavOro, #0 ) at x = L/2 , T = L, 



(III.5.3) 



2 
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0.005 0.01 0.015 



(a/L/ 

Figure 29: Continuum extrapolation of <r m = \im a /L-^o and <7 m . Two different discretizations for 
the static quark action are used and extrapolated to a common continuum limit [140]. Dividing by L2, 
the j/-axes covers about an energy range of 300 MeV. 



with (see Sect. III. 2. 3) 

f* v (x , e) = 1 in [ - f A ( Xo , e) (k v (x , e)f] . (m.5.4) 

However, at order l/m b the energy T av depends on c^ ET ,Cy^ ET . This is clearly 
inconvenient and can be avoided by choosing instead 

ri(L, e ) = - 9t + 9 t '. Fl ( L7 do ) at r = l/2 , (in.5.5) 

with 

F 1 (L,0) = iln[/ 1 (0)(A :i ^)) 3 ]. (III.5.6) 

A spin average is taken which will be relevant when we include the first order in 1/ m b ■ 
Both r av and T\ turn into m|j v when L and T are large. 

In the numerical evaluation we chose 9q = 0. We require matching, $? CD (Li,M b ) = 
$f QET (L!,M b ), with 

$? CD (L!,M b ) = 6<o) , (III.5.7) 

$f QET (L l5 M b ) = L 1 (rf at (L 1 ,c? )+m b ). (III.5.8) 
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Figure 30: Continuum extrapolation of Mb), for z = 10.4 , 12.1, 13.3 from bottom to top [140]. 

On the right, the equivalent in the alternative strategy is shown with 9o = 1/2. Dividing by Li, the 
y-ax.es covers about a range of 2 GeV. 



Here T| tat refers to eq. (III. 5. 5) at the lowest order in l/m^. Eq. (III. 5. 2) can then be 
written in the simple form, 

d>f QET (2L,M b ) = 2^ ET (L,M h ) + a m (f(L)) , (III.5.9) 

a m [g\L)) = 2L[Tf a \2L,9o)-rf a \L,e )}. (III.5.10) 

In cr m the divergent 5m as well as the mass shift m\y cancel. Its continuum extrapolation 
is illustrated in Fig. 29. 

We now see immediately that 

m ™ = mf ) +mf ) + 0(A^ CD /m b ) (III.5.11) 
mf ] = E stat -rf*\L 2 ,e ) (III.5.12) 

S v ' 

a-o in HQET 

m (0a) = r stat (L2 ^ o) _ r stat (Li ^ o)+ 1 $ QCD (Li)Mb)) (ffl .5.13) 

a-^o in HQET a-,0 in QCD 

where the first term in eq. (III. 5. 13) (in dimensionless form) is given by eq. (III. 5. 10) 
and -E s tat is the infinite volume energy of a B-meson in static approximation introduced 
earlier. It is often called the static binding energy. As indicated, the continuum limit 
can be taken in each individual step; a numerical example for the last term is shown in 
Fig. 30. 

After obtaining all pieces in eq. (III. 5. 11), the equation is numerically solved for 
Zb = MbLi, see Fig. 31. Since also the size of L\ in units of ro [47] is known, one can 
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Figure 31: Graphical solution of eq. (III. 5. 11). Data points are 2$i and the horizontal error band is 
i 2 m| v - a m - L 2 [E Bt ^ - If at (L 2 , 6 )]. 

quote 

r A*£ 0) = 17.25(20) -> Aff } = 6.806(79) GeV , (III.5.14) 

where ro = 0.5 fm is used. This result is in the quenched approximation but includes the 
lowest non-trivial order in 1/mb and a continuum limit. In the form of eq. (III. 3. 22), 
we have defined the subtraction E^ t = r| tat (L2)- The term rrig"' is then given entirely 
by the finite volume computations while mg*' results from the large volume i£ s tat an d 
r| tat (Z/2). We proceed to discuss the 1/mb corrections. 

III. 5. 3. 2 Including l/m^ corrections. 

Both the spin average m B v and the finite volume energy Ti are constructed such that 
a; sp i n drops out of their 1/mb expansion. Furthermore, Ti does not obtain any 1/mb- 
terms due to the expansion of composite fields, a convenient property of the Schrddinger 
functional boundary fields. Thus, for the particular problem of relating m|j v to the mass 
of the quark, only one more matching observable has to be defined, to determine Wkin- 
The choice in [140] was 

$ 2 (L,M b ) = R 1 (L,e 1 ,e 2 )-Rf at (L,e 1 ,e 2 ), (in.5.15) 

R 1 {L,9 1 ,e 2 ) = Fi(L,0i) -Fi(L,0 2 ) atT = L/2, (III.5.16) 
Rf at (L, 9 U 6 2 ) = In [/f 8 *^!)//? 8 *^)] , at T = L/2 (III.5.17) 

where the static piece i?| tat is subtracted such that $2 is proportional to ^kin- It is now 
a matter of simple algebra to compute the step scaling functions and finally the 1/mb- 
correction to eq. (III. 5. 14). In the numerical evaluation three combinations (61, 2 ) = 
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r M b = r (M^ 0) + M^' 








6»i = Si = 1/2 
6> 2 = 1/2 6» 2 = 1 


6»! = 1 
02 = 





17.25(20) 


Main strategy 
17.12(22) 17.12(22) 


17.12(22) 



1/2 
1 


17.05(25) 
17.01(22) 
16.78(28) 


Alternative stratej 
17.25(28) 17.23(27) 
17.23(28) 17.21(27) 
17.17(32) 17.14(30) 


17.24(27) 
17.22(28) 
17.15(30) 



Table 5: RGI results of M b in the static approximation and of the l/m b correction for 
the alternative strategy. 



(0, 1/2) , (1/2, 1) , (0, 1) were chosen. The l/m b -correction, split as in eq. (III. 3. 22), but 
for the quark mass, turned out to be about independent of 61,62 with (after continuum 
extrapolation) 

r M^ a) = -0.06(3), r M^ lb) = -0.06(8) -► r M b = 17.12(22). (III.5.18) 

Indeed since for our choice of matching condition the static approximation is indepen- 
dent of 61,62, the values for r$M^ have to be independent up to small fo ^-QCT>/ m b 
terms. This is so because the result including all terms has this precision. 

With A^jgro = 0.602(48) [17,93], the 4-loop (3 function and the mass anomalous 
dimension [174-177], we translate M b = M^ 0) + to the mass in the MS scheme, 

m b (m b ) = 4.347(48)GeV ; (HI.5.19) 

the associated perturbative uncertainty can safely be neglected. In the MS scheme the 
l/m b term amounts to — 27(22)MeV. 



III. 5. 3. 3 An alternative strategy. 

Also an alternative strategy has been tested. It is based on T av , eq. (III. 5. 3). It hence 
also involves the l/m b -corrections to the currents, but only in the combination c^ ET = 
(c^ ET + 3cy^ ET )/4. Three matching observables are necessary. They can be found 
in [140]. Here we only mention that again three angles 60,61,62 appear, but now in 
nine different combinations. Together with the above discussed strategy, there are 12 
different matching conditions. 

At lowest order in l/m b , the equations relating quark mass and the spin-averaged 
B- meson mass remain unchanged. Only the step scaling function er m is replaced by 
<7 m , and $1 is replaced by $1. Both depend on the angle 6$, but not on 61,62. A 
comparison is shown in Fig. 29, Fig. 30. Dividing by L\ and L2 respectively, one 
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observes differences of about 250 MeV in these quantities. According to the HQET 
counting they are expected to be of order Aqcd- In the combination, which leads to 
, these differences are reduced to 0(Aq CD /mb). 

In the table shown, such differences are barely visible in the lowest order results, 
where one would estimate 0(Aq CD /m^) ~ 1% effects. After including the l/m^ 
corrections no signs of differences remain. Indeed one would estimate them to be 
0((AQCD/"ib) 3 ) ~ 0.1%, quite below our statistical precision. We conclude that the 
expansion works very well, as expected. 

We finally note that the overall uncertainty of the b-quark mass is dominated by 
the renormalization factors in relativistic QCD, i.e. by the relation between the bare 
mass and the RGI mass in QCD. The second most important contribution to the errors 
is the large volume E^ in . This quantity contains a quadratic divergence. Once 
is subtracted, the remaining finite piece has errors which grow fast as one approaches 
the continuum limit. However, the methods of [178] have not yet been applied and we 
expect to reduce this source of error. 

III. 6 More literature 

Before summarizing and discussing the perspectives for HQET on the lattice, let us 
attempt to give a brief guide to literature on subjects which we did not cover. 

Early numerical computations have been summarized in reviews [162, 179, 180]. 
The power divergent mixing of operators of different dimensions was seen in an explicit 
perturbative computation [181]. The B* — B mass splitting, which is a l/m^ effect given 
at lowest order by E sp i n has been investigated with perturbative w S pin by [182-185]. 

In static approximation, the four fermion operators responsible for B — B mixing 
with Wilson fermions require the subtraction of operators with different chiralities, 
which was attempted with perturbatively estimated mixing coefficients [184,186]. Later 
it was realized that the mixing is completely avoided with an action with exact chiral 
symmetry (and spin symmetry) [187] and a calculation of B — B mixing was carried out 
with perturbative renormalization [188]. 

The dependence of heavy-light meson properties on the mass of the light quark 
can be described by a suitable chiral effective theory [189, 190]. Its Lagrangian involves 
the B*B7r coupling as a low energy constant. This was estimated in a static compu- 
tation [191]. Also an exploratory computation of the Isgur-Wise functions ti/ 2 , T3/2 at 
zero recoil has been carried out [192, 193]. These functions give the form factors for 
transitions between heavy-light mesons of different orbital angular momentum. 

Finally, the formulation of HQET at a finite velocity has been investigated by Agli- 
etti et al. We refer to [138,194] and references therein. Both perturbative investigations 
of the renormalization [193, 195] and numerical simulations [192, 196] have been carried 
out. 
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III. 7 Summary and perspectives 

Non-perturbative HQET at the leading order in l/m-b has reached a satisfactory status. 
To underline this statement, we comment briefly on the progress made in recent years. 
The ground breaking work of [132, 197, 198] was followed by intense activity leading 
to a computation where (within the quenched approximation) all error sources were 
investigated and controlled as well as possible in the middle of the nineties [199]. For 
the B decay constant, the authors estimated three sources of errors between 7% and 12% 
each, one of them due to an estimated precision of the perturbative renormalization. 

It took almost a decade until a non-perturbative method for the renormalization in 
the effective theory was fully developed [9, 10], but meanwhile a total error of 4% in the 
B s -meson decay constant has been reached [200], including the continuum extrapolation. 
In reaching this accuracy also the reduction of statistical errors [145, 146] was relevant. 

Also the standard model B-parameter for B-B mixing is well on its way [11,158,201]. 
In these cases, bare perturbation theory is now avoided by non-perturbative renormal- 
ization of the lattice operators. As explained in Sect. III. 3.1, all static-light bilinears 
require only one common renormalization factor Z^^ Gl . As this is known [10], also 
semi-leptonic decay form factors can be computed with non-perturbative renormaliza- 
tion. For completeness we note that a source of a perturbative error remains in Cps, 
Fig. 24, and its relatives for other bilinears, but here high order continuum perturbation 
theory is available [150-153] and is applied at the b-scale. This error is under reason- 
able control. Applying these methods to the theory with dynamical fermions is straight 
forward; "only" the usual problems of simulations with light quarks have to be solved. 

By themselves such lowest order (in results are not expected to have an in- 

teresting precision for phenomenological applications, but certainly they can constrain 
the large mass behavior computed with other methods [161-173]. An interesting ap- 
plication has been the combination of the approach of [172, 173] with the large mass 
behavior computed via HQET [200]. 

The 1/mb corrections can also be computed directly in the effective theory. Here, 
the necessary steps have been carried out in detail for the mass of the b-quark in 
the quenched approximation. The resulting precision is rather satisfactory and higher 
order 1/mb-corrections can be neglected. The latter has also been verified explicitly 
by comparing the results following from a number of different matching conditions. No 
significant differences were found. 

An extension of this computation to full QCD has been started by the ALPHA 
collaboration. This is of particular interest because for full QCD it is of course more 
difficult to reach the small lattice spacings needed for computations with relativistic 
quarks of masses around m c h a rm ~~ the basis of the alternatives mentioned above. It 
is time to address other observables such as the B-meson decay constant in this direct 
approach! 
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